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ERRATA 


Automation and Remote Control, Vol. 22, No. 1, translation published August, 1961 

On p. 46, 2nd paragraph, 2nd line, instead of “any values of A; and Bj “. it should read “any values of A, and 
By." 

On p. 87, 4th paragraph (1st paragraph extending full width of page), next to last line, instead of “conditions 
(26)," it should read “conditions (28). 


REPRESENTING NONSTATIONARY LINEAR DIFFERENTIAL 
POLYNOMIAL OPERATORS IN THE FORM OF SUMS 
OF STATIONARY OPERATORS 


A. N. Sklyarevich 


(Riga) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 3, 
pp. 297-305, March, 1961 

Original article submitted October 1, 1960 


It is being shown that an arbitrary nonstationary linear differential polynomial operator can be repre- 
sented in the form of a sum of stationary operators with associated kernels, An algorithm is given which 
enables one to find all the components of the representation. In the Appendix a formula for transform- 
ing of linear operators is derived of which use is made in the theoretical investigations and in the de- 
termination of the structure of automatic control systems, 


In many cases the study of the structure of an automatic control system, the determination of the weighting 
function as well as establishing the exactness criteria of the operation of the system can be simplified considerably if 
its behavior, describable by an equation 


O (D, t)z = f(t), (1) 
with & (p. t) being the known nonstationary operator, x(t) the sought response to the excitation f(t), is transformed into 


> ®, (Y) (9, @ =) 
k=1 (2) 


In the latter equation & (p) are the known stationary operators, and j (t) are known functions of time. 


The case of s = 1 is a particular one of the formula (2). Now the study of the nonstationary system can be re- 
duced to that of a stationary one, It is certainly of interest to show that such a reduction is actually possible, and to 
provide an algorithm which brings the equation (1) to assume the form (2), When s >1 the weighting function and 
the exactness criteria may, in a number of cases, be obtained in a fairly straightforward manner, if use is made of the 
Laplace transforms. 


The operator & (p. t) will be regarded as differentiable if it can be expanded into a Taylor series with respect 
to the differentiation operator p= d/dt. The operator & (p. t) will be regarded as nonstationary if it depends explic- 
itly on the independent variable of differentiation, and will be regarded as stationary when there is no such depend- 


ence, We shall presume that a nonstationary operator & (p. t) is reducible to the stationary & (p) with kernel ¢ (t) if 
for any sufficiently differentiable function x (t) there is valid the equality 


(v, t)z = ®, (p) (9 (2) 2). 
The reducibility of the operator & (p. t) to a stationary operator & (ꝓ) with kernel (t) will be denoted by 
(v. t) = O(p) f). 
. e., cutting out, for brevity, the arbitrary function x(t), as well as the closing bracket. 


We shall show that an arbitrary nonstationary linear differential polynomial operator & (p. t) which can be ex- 
panded into a finite Taylor series with respect to the operator p and the independent variable t can be represented in 
the form of a sum of stationary operators with associated kernels: 


(. t)= Lo. (P) (t). (4) 


(3) 


— 

— 

— 


We shall establish as a preliminary result the conditions under which an operator & (p.t) reduces to a stationary 
operator ©,(p) with kernel ¢ (t), These conditions can be found by using the formula 


®, (v) l (t) = p(t +9) (p) = Op, 5), * 


in which qs 0/0p , and whose derivation is given in the Appendix. A particular form of the formula (5) without der- 
ivation and without shortened notation was given in A. V. Solodov's work [I]. 
We expand the operator v (t+q) as a Taylor series with respect to the operator q: 


(r) 
r=0 (6) 


and substitute this result in (5) and use the identity 
7 ®, (p) = Of" (p). (1) 
As a result, we obtain 


2 % 
O(p, t) = +9) — 2 
— (8) 


An arbitrary nonstationary linear differential operator & (p. t) reduces therefore to a stationary operator provided 
there exist functions v (t) and & (p) which satisfy the equation (8). 


We now determine the functions ¢ (t) and & (p) assuming that such a reduction is possible and also that & (p..) 
is a polynomial of degree m with respect to the variable t and of degree n with respect to the operator p: 


t) = Dd) cistip’. 
imo j=0 


(9) 
In this case it seems reasonable to take 
k=0 l=o (10) 
We represent the functions ¢ (t) and & (p) in the following form: 
m n 
1 
k=0 
and we designate 
Ry = min Im — k. I]. ay 
Then with r Ny we have 
m n 
(t — kl a 11 b,p'—; 
) 23 1 (Pp) ip 
but when r > Rog, one of the derivatives of order r vanishes identically, Consequently 
Re m un 
i—r 
We introduce new indices i = k—r, j = I- and invert the order of summation in the formula (12): 
min Rij 
i 10 
0 222 (13) 


i=0 j=0 r=0 


— 


(5) 


(6) 


», t) 


(9) 


(10) 


(11) 


(12) 


(13) 


Comparing (13) with the equality (9) we find 


Rij 
r=0 


When the coefficients cjj are known, the relations (14) represent equations with unknown products ajbj. We in- 
troduce the notation 


= dy (15) 


and it can be shown by the method of mathematical induction that the solutions of the equations (14) can be expressed 
by the formulas 


Rij 
d= > Cir, 
r=0 


(16) 


With the known matrix ||cjjll of the coefficients, the equations (16) determine the matrix of the quantities di il. 
Conversely, when the matrix ¶dijll is known, the formulas 


r=0 


(17) 
provide the corresponding coefficient matrix cijll. The matrices leijll and Hdijll have an equal number of rows and 
an equal number of columns, It proves helpful to number their rows from bottom upwards and the columns from left 


to right. The elements will be correspondingly equal to one another (in particular mn = dmn * O) in the top rows 
and in the marginal right columns of the matrices ||cjjl|_and idijll. 


Thus we can establish the (m+ 1)(n+ 1) relations 
a,b; — di;, (18) 


connecting the m+ 1 coefficients of the kernel & (t) and the n+1 coefficients of the operator 1 (p). When the oper- 
ator & (p. t) is reducible to the operator & (p) with kernel ¢ (t), all the above relations are true. Taking the particu- 
lar values i=a, i=8, j=y, j=5, we obtain 


a,b, = dey, = das, agb. den, agbs = dps, 
and hence 


As the indices a, 8, y, 5 can assume any integral values within the bounds OS m, 0=8=m, 0= y=n, 
oo <n, the relation (19) is valid for all the second order minors of the matrix Idijll. This signifies that the ma- 
trix Nj jl is of rank one [2]. 


Therefore, when the operator & (p. t) is reducible to the operator o (p) with kemel ¢ (t), the matrix I diſll. cor- 
responding to the coefficients of the operator & (p. t),is of rank one, The converse statement is also true: if the ma- 
trix I dijll corresponding to an operator & (p. t) is of rank one then the operator & (p.t) is reducible to an operator 


(p) with kernel (t). The proof of the converse statement follows from the simultaneous validity of the equations 
(18). 


We determine the coefficients aj and bj in the case when the operator is reducible. The representation 
(p, t) = (p) (e) 


is obviously not a unique one, as all the coefficients aj can be increased k times provided that at the same time we 
reduce the bj coefficients by an equal number of times. One of the coefficients aj or bj can, therefore, assume an 
arbitrary value, It seems in order to take bp = 1 in view of the usual way of listing transfer functions; then 


din = abn = d (20) 


— 
1 
| | 
— ̃ͤ 
— 

=| 


We shall now concentrate on the value i= m. As dmj = ambj. 


4 
(21) 


But Cin = din · emj = dmj: therefore 


(22) 


The formulas (20)-(22) determine the coefficients of the operator & (p) and of the kernel & (t) when the oper- 
ator & (p. t) is reducible to the operator S1 (p) with kernel & (t). All these coefficients only depend on the quantities 
in the top row and in the furthest right column of the matrices leijll and I dijll. 


It follows from the above that in the general case a nonstationary operator & (p. t) is not reducible to a station- 
ary one ,(p). If an arbitrary nonstationary polynomial operator is given, & (p. t), then by using only the top row and 
the last right column of the coefficient matrix ||cjjll_of the operator & (p. t) one is able to determine the correspond 
ing operator ,(p), the kernel Hit) and the matrices llcjjqll , Idijill of the operator & (p) lt). In this case, how- 
ever, the identity & (p. t) (p) ft) is not valid. We shall now find the operator & (p. t) = & (p. t)— ©, (p)[¥,(t). 


All the elements vanish in the top row and the last right column of the coefficient matrix of the operator 
,(p,t). This means that the operator 1 (p. t) represents a polynomial of the m-1 degree with respect to the vari- 
able t, and of the n-1 degree with respect to the operator p. Applying this procedure Rog~ 1 times we arrive at anon- 
stationary operator vm (t) if msn, or if m >n at a stationary one Sn (p). In the coefficient matrices of the operators 
m (t) and n (p) the nonvanishing elements can only be found in the bottom row and the first left column. But, as 
the nonstationary operator v m (t) reduces to the stationary operator pꝰ with kernel Vm (t). therefore by a successive 
adding up of the relations 


Dy (Y, t) = Dea (Pp, f) — Or (v (4), 


which express the outcome of each transformation stage, we obtain the formula 


(23) 


Re 
o (v. t) = a ®, (p) (9, (t)- (24) 


The above is a formula which shows that it is possible to reduce an arbitrary nonstationary polynomial operator 
to a sum of stationary operators with associated kernels, 


In view of the linearity of the formulas (16) and (17) with respect to the quantities dij and cj;, one can give 
the following simple procedure to determine the operators p) and the kernels (t). The matrix |ldjjl]_ is found 
from the coefficient matrix ||c;;l] of the operator & (p. t). The former matrix is represented in the ſorm of a sum of 
matrices of rank one in such a way that by adding the first k matrices we obtain a matrix with k upper rows and k 
columns furthest on the right, which are the same as in the matrix dijll. From the k-th matrix we detach the top 
row and the last right column both of which consist of nonvanishing elements, Taking into account that this matrix 
represents the matrix I di jxll corresponding to the operator & (p) ſ oH (t). the coefficients bj, of the operator 4, (p) 
and the coefficients aj, of the kernel Y (t) are found with the aid of the formulas which differ from (20) and (21) 
only by their index numbering. 


Examples 
1, The behavior of a dynamical system is described by the equation 
( 4 + 2 (24 + + 6) + ( + 65 + 7) =f (t), (25) 
where f(t) is a known function, In this case 
t) = (1 + 2 + + (6 + + 20% p+ (7+ E+ (26) 
and consequently m = n = 2, 


1. 


(21) 


(24) 


The coefficient matrix ||c;jl| of the operator & (p. t) is now 


1214 
Neg l= 6 8 2 
761 
Using the formulas (16) we find the matrix dijll: 
12 1 | 
124 


The matrix N¶dijll is of rank one. Therefore the operator & (p. t) reduces to a stationary operator (p) withthe 
kernel ¢ (t), and in agreement with equalities (22) we have 


®; (p) = 1 + 2p + N. 
(t) = 1 + 2t + 2. 


The equation (25) is thus equivalent to the equation 


(p® + 2p + 1) % + 6% = / (0). (27) 
2. Given a nonstationary operator 
D(p, 1) = — 34%) p* + (84 4 31 — 10¢ — 1) p + (4t® + 124% — 2t — 10). 


It is required to reduce it to a sum of stationary operators with associated kernels, Writing down the coefficient ma- 
trix lleijll of the operator & (p. t): 


4 6 2 

210 0 
—10 -1 0 


we find the matrix dijll with the aid of the formulas (16): 


4 6 2 
4 2 0 
0 
The matrix dijl is of rank three which shows that the operator & (p. t) cannot be reduced to a single station- 


ary operator ,(p) with the kernel ¢(t), The matrix dijll shall be presented in the form of a sum of matrices of 
tank one such that the rule on obtaining the elements of the original dijll matrix is satisfied: 


4 6 2 4 6 2 0 00 00 0 
00 0 
4 2 0 0 0 0 4 20 00 0 
0 0 0 —2 —1 0 — 0 0 


To each of the component matrices there corresponds a definite stationary operator d (y) and also a definite 
kernel (t). Taking into account the formulas (20) and (21), we find 


O(pP)=2+3p+ = 2° — 
= 2+ p, G2 (t) = 21 —4, 
®;(p) = 1, Ps (t) = —4. 


= 

(22) 
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2 

we 

(t). 

ri- 

a non- 

rators 

as 

ive 

(23) 
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Consequently, for any function x(t) we have 
O(p, r = (p* + 3p + 2) [(2t® — 30 x} + (p + 2) ((2t — 1) 2] — 4 


It would be of interest to mention that in the investigated cases the rank of matrix [ld;;||_ is equal to the num- 
ber of the components in the representation of the given nonstationary polynomial operator by a sum of stationary 
operators with associated kernels, 


APPENDIX 


Transformation Formula for Linear Operators 
We shall derive a formula which will enable one to simplify the performing of transformations of the type 
y= O(p, (28) 
where (t) and x(t) are functions of the independent variable is the former assumed known, the latter unknown and 
sought. 
We shall assume that there exists a Laurent expansion for the function ¢ (t) in powers of the variable t. and also 
for the operator & (p. t) in powers of the operator p: 


= Yat, 
i=—oo (29) 


@(p, 25,0 


—00 (30) 


(, i) >} (. (2). 


We now find the value of the component & (p. t) (tix) for different values of i, Wheni=1 and j is an integer, 


we have 
p? tz) = tpiz + ip 17. 
(tx) U ( 


We multiply the- equality (32) by the function bj (t) and then we sum with respect to index i. From the formula 
(30) the relation 


a® (v. t) 
> jb; (t) op 
j=—oo (33) 


follows and the result of the summation can thus be written as 


O (p, t) (tz) = (p, t) + —— *. (34) 


The formula (34) is a generalization of a formula known in the operational calculus [3]: 
(p) (tz) = Ito (p) + (p)] z. 


By introducing the differentiation operator q 9/8p we are able to write the formula (34) as 


O (p. t) (tz) = [(t + la. (35) 


By twice applying the formula (35). we find 
(. t) (t) = It + g) O (, (tz) = (t+ 9) 9) 


or more briefly, 


(Y. i) (h = [(t + O(p, 


— 
Also 
12 (31) 
28 
260 


Analogously, we establish 


0 (p, i) (ltr) = [(t . 


(36) 
where i denotes a positive integer, 


9 The equality (36) is also valid when 1 = 0, simply from the way it is written down. This equality can be trans- 
formed by raising formally the sum t: q to the i-th power with the subsequent use of the formula 
(37) 
(28) 
and 
also 
(29) 
z 1 
(30) 
A generalization of the formula (38) is 
= i 
which is proved for j O by mathematical induction. 

* We assume that the formula (39) is valid when j=n. Then 
(32) 
rmula 

(33) 4 

3⁴ 1 23 —1 

After some reduction and by introducing the operator q we obtain 


— 
) 
) 
or 
261 


Thus the formula (39) is valid when j = n+1, and consequently it is valid for any j O. A direct check will 
show the correctness of the formula when j=0. We shall now study the values of j< 0, By applying successively the 
integration by parts formula 


1 1 
we obtain, when j=~-1, 
(an 


If the function x/t can be expanded into a Taylor series, the last term on the right-hand side of the formula 
(41) becomes arbitrarily small with n=. Therefore, when proceeding to the limit, we can suppose that 


+...) 


A generalization of the formula (42) is 


aa 


which can be proved by mathematical induction, The representation of the operator Es : 1 I in the form of a 
series with respect to the operator P. the subsequent application of the formula (40) to each term of this series, the 
introduction of the operator q and the contracting of the series to the equivalent operator U t+? 0 | are the 
main stages of the proof, 


The formula (43) can be regarded as a particular case of the formula (39) if the latter is widened to negative 
integral values of the index I. 


We multiply the relation (39) by the function bj(t) and then we sum with respect to the index i Shifting the 
place of the function bj(t) and of the operator 1/(t+q) and taking into account the formula (30) we find 


© (p, = (44) 
The formula (44) leads to the formula 
(. i) = It + Y. (45) 
valid for all negative i, The right-hand side of the formula (45) is convergent if the function tix has a Taylor ex- 


pansion, 


The formula (36) may be looked upon as a particular case of the formula (45) if the latter is widened to in- 
clude all integral values of the index i, Transforming the relation (31) with the aid of the formula (45), we obtain 


—00 


In accordance with the formula (29), we have 
E 
alt+ =e +o) 


i=—oo 


(42) 
;- 
p 
| — 
| 
262 


and therefore 


D(p, i) (p(t) 2) = [p(t +4) O(p, z. 
By using another notation, the formula (46) becomes 
(40) 
® (p, t) (= 0 + (y. (47) 
and this shows distinctly that the component ¢ (t+ q) represents an operator acting on the operator & (p. t). In order to 
dispense with the square brackets, it is necessary to expand formally the operator ¢ (t+ q) in a series of the powers of 
(41) q, and then use the equality (37). 
Bs We should like to mention two results of the formula (46). 
1, The shift theorem 
ꝙ (p. t) = (p, =e! @ (p + a, 2. 
2. The interchange of the operators v (t+ q) and & (p. t): 
= ),. 2. 
42 
The formula (5) is a particular case of the formula (46). 
A. V. Solodov, Elements of the Theory of Linear Systems with Variable Parameters [in Russian] (Izd. AIA im. 
(49) Dzerzhinskii, 1958). 
2. V. I. Smirnov, Course in Higher Mathematics [in Russian] (Gostekhizdat, 1951) vol. 3, pt. I. 
if a V. A. Ditkin and P. A, Kuznetsov, Handbook of Operational Calculus [in Russi an] (Gostekhizdat, 1951). 
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tive ON THE SYNTHESIS OF LINEAR AUTOMATIC CONTROL SYSTEMS 
D. I. Gladkov 
the 
(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No, 3, 
pp. 306-313, March, 1961 
(44) Original article submitted April 20, 1960 
A method for synthesizing linear dynamic systems with constant and variable parameters is considered. 
(45) The structure and the parameters of the compensating networks are determined. A method for an ac- 
ex- curate realization of the integrating factor is proposed, 
The block diagram of an automatic control system is basically determined by the purpose of the system and the 
12 system's theoretical ability to solve the assigned problem. The assignment of any additional properties to the system 
) 


is realized by introducing compensating networks. For this, the synthesis problem can be stated in the following man- 
ner, 


Assume that the linear dynamic system shown in Fig. 1a is given. 


It is necessary to determine the block diagram and the parameters of compensating networks which would bring 
the dynamic properties of the given system closer to those of a certain optimum system, The dynamic properties of 
the optimum system can be conveniently characterized by a weighting function, For linear systems that are described 


- 


by differential equations, this function can be represented in the following form: 


be (t. 1) = 
j=1 


The optimum weighting function can be determined on the basis of different principles [I]. 
Let us introduce the compensating networks in the manner shown in Fig. Ib. It should be noted that 


t 
8, (f. r) = den (n. v) gt, n) dn. 


The differential equations describing the — of systems with the weighting functions gu (t. 1) and gult. j 
will be given by 


tou = Zin, 

pa (2) 
>) = Zins 
— (3) 


respectively, where p = d/dt, the exponents p and p determine the order of the derivatives, and k and m determine 
the order of the differential equations, The coefficients ap and by are generally functions of time. 


Let us determine the weighting function gu (t. 1) of the compeusating networks, 


Let gf (t. 1) be the weighting function of the forward network in the block diagram of the assigned automatic 
control system with the introduced compensating networks, and let g (t. ) be the weighting function of a series as- 
sembly of compensating networks and of a dynamic system with the weighting function gi (t. T). 


For determining gf (t. 1), it is necessary [2] to solve the Volterra integral equation of the second king 


Bx (t, = Bolt, 0) + gi b. Je. (Gs an 


* 


(4) 


Equation (4) can be solved by using the method of successive approximations, The get. 1) and g (n. &) func- 
tions are represented in the form of a sum. Each term of the sum is equal to the product of two factors, one of which 
is a function of the first, and the other is a function of the second argument of the weighting function. This fact en- 
ables us to state that the solution of Eq, (4) is given by 


= c(t) d. (0). 
121 


g, (. v) 


The value of n is determined by the required accuracy with which Fd. (4) is to be solved. 


— — 
a 
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If we know the gf (t. 1) function, we can determine gy; (t. 1). For this, the method of inverse factors [2] must 
be used: 


t 
eu (t,t) den (n. dn. ©) 


By using the known gui (t. 1) and gis (t. 1) functions, the weighting function of the compensating networks is de- 
termined in exactly the same manner: 


t 

6. 9 = Ven +) g(t, db. 

In Eqs. (6) and (7), gũ (t. 1) and ga (t. ) denote the weighting functions of systems that are inverse with respect 

to the dynamic systems described by differential equations (2) and (3), respectively, 


The functions gũ (t. 1) and gi9(t,T) can be determined by putting xout = 6(t—T) in equations (2) and (3), With 
this, we obtain 


By (t. t) = Nas (i) det — zy, 
* 

en (t. t) = Dj b, (t) (¢ — x), 
»=0 


By substituting (5) and (8) in Eq. (6), we find 


tm n 
* p=0 1 


n mnt 
121 92117 


Here, the notation di (r) = a9(t) dj (rT) is used. 


By further substituting (9) and (10) in Ed. (7), we obtain the expression for the weighting function of the com- 
pensating networks: 


n 


da (x) + 


11 


mn * 
7 p=0 


f which or 
act en- n mnt 
Bx t. t) = Tos (t) de (t) + 2 2 fain — T) as (n) Cot (t) di (n) dn + 


ke 
dos (x) dy + an 


denn — ay (n) (E) di (n) (t) &) dn dt. 


(3) 
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(9) 
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Here, the notation i (t) = belt) ei (ti) was introduced. 


An analysis of expression (11) indicates that the compensating network represents a parallel combination of dy- 
namic systems with weighting functions of the following types: 


en (ts = cat dn 
Boa (ts = die (en (0) de an, 


t 
Boo (t, Je. dos (x) by. a, 


deen n ay () (E) dn) (t) — &) dy dt. 


Boa (. 


We shall now see how systems with such weighting functions can be realized. 


The block diagram of a dynamic system with the weighting function gu (t. T) will be composed in the form 
shown in Fig. 2a [3]. 


The integrating element is preceded by the amplifying element with the variable gain dgj (t). which is equal to 
that portion of the weighting function which depends only on the second argument; the integrating element is fol- 
lowed by the amplifying element with the variable gain cg (t). which is equal to that part of the weighting function 


gu (t. 1) which depends only on the first argument. 
An elementary check: 


t 
de (i () (nan = de (2), 
doi (x) co (t) = Bor (t, t) 


confirms the above reasoning. 
Let us compose the block diagram and determine the parameters of a dynamic system with the weighting func- 


tion gog{t,T). 


2E LHA 


0 d 


Fig. 2. 


On the basis of the expression for the weighting function pertaining to the series coupling of a system or its ele 
ments, the following can be stated. A dynamic * with the weighting function gi (t. 7) represents a series as- 
sembly of systems with the weighting functions 

Boas (t, 2) = \ den Ya, and 


tt 


gon (t, T) = by (t) 5 (t — 5). 


＋˙ũ„“5 — 
11 — 
— —V 
11 


ts ele- 


A system with the weighting function g (t. 71) is represented by a series assembly of systems with the weight- 
ing functions 


t 
(t, t) = jo n — T) (n) cy (t) di (n) and Boss (t, c) 5) t). 


A system with the weighting function geit, T) represents, in turn, a series assembly of dynamic systems with 
the weighting functions 


Boars (f, t) = det — c) and gi (t, c) = ei (t) a, (c der. 


Figure 2b shows the block diagram of a dynamic system with the weighting function gas (t. 1). 


The block diagrams of systems with the weighting functions ge (t, 1) and gu (t. 1) can be obtained in a similar 
manner (Fig.2,c and d). 


Thus, a compensating network which would bring the dynamic properties of the given system closer to the opti- 
mum properties will represent a parallel assembly of networks consisting of differentiating, integrating, and amplify- 
ing elements with variable gain values, 


A different method of composing the block diagram of compensating networks can be used, 


In analyzing Eqs. (6) and (7), the following can readily be noticed, A system with the weighting function 
8k (t. 1) represents a series assembly of systems with the weighting functions gj; (t. 1) and gi (t. 7). On the other hand, 
a compensating network with gu (t. 1) represents a series assembly of systems with the weighting functions g (t. 7) 
and gi(t,T). Consequently, the compensating network consists of a series assembly of systems or elements with the 
weighting functions gi (t. 1). gf (t. 1). and gig(t,T). Hence, by the way, the following rule: The weighting function 
of the middle of the three elements connected in series (see the forward path in Fig. 1b) is determined by substituting 
the weighting function gf(t,r) of the entire assembly for the g (t. 1) function which is to be found, and by replacing 
the weighting functions of the end elements by the corresponding functions of the inverse elements, 


7 * By taking into account Eqs. (5), (8), and (9), we obtain the block dia- 

82 te 5 3 gram of the compensating network (Fig. 3). 

8 t In constructing actual compensating networks, the realization of the 
A 1 riuft integrating element entails considerable difficulties, A method of accurate 

the i 2 € replacement of the integrating element by an inertial element with vari- 

b able parameters [3] is proposed below. 

11 a(t) 1 700 Io The weighting function of an inertial element with the time constant T 
pT is given by 
c 


g(t, re er. 


ual to 
* 
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The block diagram of such an element can be represented in the form shown in Fig. 4a. If we connect ampli- 
fying elements whose variable gains are respectively equal to re UT and e/T to the input and output of such an ele- 
ment, we shall obtain an integrating element (Fig. 4b). 


The block diagram of a system with the weighting function go; (t. 1) can be now represented in the form shown 
in Fig. 4c, where the following notation has been introduced: 


gelt) = (t) e Tt, = cos (t) ed. (19 


In practice, such a system can be realized in the form of an RC-filter (Fig. 5) with a variable resistor and a 
variable capacitor, 


The resistance R(t) is determined according to the equation 


doi 


He 


The capacitor variable capacitance is given by the equation 


t 
: an 
(14) 


Equations (13) and (14) are derived on the basis of the following principles, By using the diagrams in Figs, 2a 
and 5, we can readily obtain the differential equations describing the operation of the corresponding dynamic system: 


— = Ci (t) dos (t) 2 
6% (4) in’ (15) 


. 1 1 1 


The following conditions must be fulfilled if the equation uout = xout is to be satisfied: 


Cos (e) 


1 


By solving the system of equations (17) with respect to R(t) and C(t), we obtain Eqs, (13) and (14). 


The ci (t) and do; (t) functions are components of the weighting function gq (t. 1). For the element in questi, 
Coi (t) = doi (t) = 1. By using Eqs. (13) and (14), we find: R(t) = et, C(t) = et. 


The considered synthesis method is applicable to steady-state as well as to unsteady-state automatic control 
systems, 


Example, The following signal acts on the output of the dynamic system whose block diagram is given in 
1 
Z()=S() 


where S(t) is the signal useful portion, and X(t) is the noise, 


The useful signal is a linear function of time, S(t) = Ut, where U is a random quantity whose mathematical & 
pectation is equal to zero and whose dispersion is Du · 


We shall assume that the interference is given in the form of white noise with the autocorrelation function 
K. (t. r) = 4% (t — 8). 


: t t 
At / t 
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The forward path of the system's block diagram contains an inertial element, the inverse weighting function of 
which is given by 


(18) 


Feedback is provided by an integrating element with the weighting function 


62 (t. t)=1(¢—T). (19) 


It is necessary to determine the block diagram and the parameters of compensating networks which would se- 


(12) 


ia cure the separation of the useful signal derivative from noise with the minimum root-mean-square error. 
The weighting function of a dynamic system which would be optimum in the above sense can be determined by 
using the well-known methods [1]. In this case, 
(13) 
It is obvious that 
Bo(t, c) Jo (t) @o(t), 
(14) 
Je ( = * 
11221 
igs. 2a + 3k, (22) 
(23) 
(15) We shall assume that the compensating networks with the weighting function gi (t. 1) are connected to the sys- 
tem in the manner shown in Fig. 1b. 
(16) In the example under consideration, gu (t. 1) = gi t. 1): gut. 1) is the weighting function of the amplifying ele- 


ment whose gain is equal to unity, It is obvious that 
Big (t. c) = b(t — 5). 

(24) 

From Ed. (4), we find 


t 


* * 


We shall assume that the first approximation provides a sufficiently 
accurate solution of the integral equation (21), Then, 


By (6. 9) = fa lt) + + 
+/s(t) Ps (t). 


et). = e), 
(t) = 7e (t) e (i). 
(T) = [Go (c) a (c) e — + 
+ % ()]. 
G2 (Tt) = . = a(t), 


„ =\ 


t 
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If we know (18) and (26), by using Eq. (6), we find 


t 
T ; 


* 


By using Ed. (7) and by taking into account Ed. (24), we find that gi (t. 1) = gt. 7). 
In correspondence with Figs, 2a and 2b, we shall obtain the block diagram of the compensating networks (Fig. 6a). 
The structural transformations are given in Fig. 6b. 


The same diagram will be obtained if the compensating networks are realized directly as a series assembly of 
systems with the weighting functions gi, (t. 7). gf (t. 1). gut. 7). 


The actual expressions for the Fi (t), (t). (t). (t). f(t), (t) functions depend on the and (t) fune- 
tions, as can be seen from the equations, In determining a(t), o (t). and c(t), the integration can be performed ana- 
lytically, numerically, or graphically, No substantial difficulties will be encountered in this. 


In the example under consideration, the a(t) and c(t) functions can be determined analytically in a relatively 
simple manner: They are given by the following expressions: 


a(t)= — + are tg 2 


1 3 


The b (t) function can be more easily determined by using the numerical or the graphical method. 
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A SERVOSYSTEM WITH LOGICAL CONTROL 


E. Kk. Shigin 


(Voronezh) 
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pp. 314-321, March, 1961 
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The present article describes an electromechanical servosystem where the parameters of the compen- 
sating elements change in steps, which is secured by means of a semiconductor logical unit. The logi- 
cal laws of the parameter jump- like variations are assigned in such a manner that a system with second- 
order astatism with respect to the control input and with first order astatism with respect to the disturb- 
ance is not subject to overcontrol for jump-like control inputs, The experimental results are given. 


There are definite limitations imposed on the performance figure of automatic control systems containing com- 
pensating elements with nonvariable parameters, For a given system to be controlled and a given final element, the 
performance figure cannot be improved beyond a certain limiting value [3], The effect of the element characteristic 


nonlinearities and the limitation of the maximum allowable acceleration of the quantity to be controlled lower the 
limiting performance figures, 


By securing jump-like or continuous variation of the compensating element parameters during the transient 
process, the above-mentioned limiting performance figures can be considerably exceeded, Certain possible laws ac- 


cording to which the parameters of the compensating elements (in the first place, differentiating elements) can vary 
have been described in the literature [2]. 


A typical block diagram of servosystems which are widely used in technology is shown in Fig. 1 [3 and 4}, A 
system constructed on the basis of this diagram exhibits second-order astatism with respect to the control input and 
first-order astatism with respect to the disturbance. 


The system consists of elements and factors (Fig. 1) which are described by the following equations, respective- 
ly: sensing element 1: 


= Ly — (1) 


differentiating compensating element 2: 


(Tip + = KyT (2) 


integrating compensating factor 3; 


or a passive integrating element: 
(Top + = Keri, 
differentiating element of the second order 4: 


(Typ +1) = K, 


adding amplifier 5: 


(T. + 1) = K. (a, + + + %), 
power amplifier 6: 


(7'sp + ). = K, 5, 
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system to be controlled 7: 


(TT. * + Top) 2, . + f (ö. (8) 
where f(t) is the disturbance, x9 is the command input, x, is the error, xz. X3, Xg, Xs, and xs are the output quantities, 
Ty. Tg, Ts, T., Ts, Tg, and T. are the time constants, and Ky, Kg, Ks, Kg, and Kg are the gains of the respective ele- 
ments and factors in the system, x; is the controlled variable, and p is the differentiation symbol. 


ft) 


Fig. 1. Servosystem block diagram. 


As was shown in (5), jump-like variation of the parameters of differentiating element 2 and integrating element 
3 must be secured in such a system in order to reduce the transient process time. 


If integration with the opposite sign is introduced together with forced damping, the transient process will take 
place without excessive overcontrol [6]. 


By securing the variability of the sign of the element 4 output quantity and the jump-like variation of Kg, we 
can reduce the transient process time and the maximum deviation of the controlled variable during the stabilization 
period [7]. 


The logical laws of the jump-like switching of the parameters of elements and factors 2, 3, and 4 can be writ 
ten in the following manner: 


(9) 


for 1143 <0 and 232, > 0, 
io for 2123 > 0, 
Ky, for 11422 and 122 C0, 
K. = K oo for 114 > 0, 
oo for C 0 and 2321 >0 


| Kx for 114 <0 and 122. C0, 
K, = 


while the equation of factor 3 will be written thus: 
K K. 


where xz is the factor 3 output value at the moment of parameter switching, 


Ky) for 2122 0 and 1224 > 0, (15) 
K., for 2,2,>0 and 22. <0, (16) 
Ky for 2122 <0. (17) 


Here, Kyp, Kg, and Kg are the basic gain values of the corresponding elements and factors. KA. Kaz, KAg, and 
Kas are the gains of the corresponding parameters, and 


Tip * = Te (K. = const, K, const). 
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Logical operations following the laws given by (9)-(17) were realized in an experimental servosystem by means 
of a unit with logical elements, to the input of which the quantities xi. xz, and xj were fed. By performing the logi- 
cal operations of comparing the signs of these quantities, the unit exerts actions A;. Az, and A; on elements 2, 3, and 
4 (Fig. 1) by varying the parameters of the latter. 


A DRK-627 motor served as the servosystem final element, the summation of signals was performed by means 
of a magnetic amplifier 5 (Fig. 1), and an electronic amplifier 6 was used for matching the load. 


Figure 2 shows the wiring diagram of the logical element unit, which is based on semiconductors. In order to 


secure the resistance decoupling of circuits in the system, the de voltages to be analyzed are transformed into pulses 
in this unit. 


The de current signals from the element 1 (potentiometric bridge), the passive differentiating element 2, and 
the second-order differentiating element 4 (Fig. 1) are fed to the inputs of the logical element unit (uini = xf. 
uin: xz, Uins x). These de signals are transformed into pulsed voltage which is in step or in phase opposition to the 
de voltage us feeding the unit in the following manner: pulses of rectified ac voltage are fed from the windings ng, ng, 
and ny of transformer Tri to the collectors of triodes TI and Tz, Ts and Tg, Ty and Tyg (selective amplifiers) through 
diodes Dg, Ds, Dyz, Dyg, Dyg, and D. For ujny = Uing = Uins = 0, pulses of pulsating voltage uz will appear at the out- 
put of the selective amplifiers (at the points a and b). as is shown in the time diagrams (Fig. 3a). The magnitude of 
the pulse peaks is determined by selecting the mixing resistors Rg, Rg, Rig. Rey, Rao, and Reg. 


The input voltages (uin i ) amplify those pulses whose 
polarity coincides with the polarity of uin i. and they suppress 
2.0 DAA 29 0 HME those pulses whose polarity is opposite to that of uin i. In this, 


20 NN. pulses ui with the duty ratio 0.5 will appear at the output of the 
0 P NN selective ampliſiers (Fig. 3a). 
n a a N Pulses u, are fed to the inputs of the relay amplifiers, which 
are based on positive-feedback two-stage amplifier circuits with 
4300 nuiodes Ts, Ty, Ty, Ts, Tu- and Tyg. The circuits and the charac 
l J \ j U a teristics of such amplifiers have been described in [1]. In contrast 


4570 “SEP to the circuit described in [1], the blocking of the output triode 
Tg(Tg, T) is secured by supplying a positive instead of a negative 


9 U 0 U 8 voltage to the relay input triode base. Such action of the positive 
* voltage is possible due to the action of the biasing resistance 
RIO (Rez, Red), which prepares the complete unblocking of triode 
T3(Tz, Ty). The collectors of triodes T. (Tg, Ty) are supplied with 


1,470 N pulsed voltage uz from windings ny. nig. and ng of the pulse trans- 


former Tr2 through diodes Dyp, and Deg (Fig. 3b). 
* f The peak value of voltage u is chosen in such a manner 


that triodes Ty, Tg, and Ty, are completely blocked for uin i C. 
For uin i 0, the suppressed series of pulses uz will not block the 
triode T. (Tg, T), and voltage uz will be applied in steps to the 
primary winding nz n. gy) of the output transformer Tr3(Tr4, Tr5), as a result of which a voltage of complex form 
will arise at the transformer Tr3(Tr4, Tr5) output, as is shown in Fig. 3c. 


From the output windings nyg and nyg, voltages up are fed to the bases of the comparison triodes Tn and Tyg. If 
the voltage at the base of triode TH or Ty, is equal to zero, the triode is blocked. If voltage up + 0 is fed to the 
base, the triode is open during that portion of the up period when up < 0, and it is even more thoroughly blocked for 
Up > 0. 


Pulses uc are fed to the collectors of triodes TA and Ty, through diodes Di and D (Fig. 3d), In the case where 
the pulses uc and the negative portions of voltage up coincide in time (the triode is open), a negative voltage is fed 
to the base of the amplifying triode Ty(Ty), whereby the latter is unblocked. The collectors of the amplifying triodes 
Tu and Ty are supplied with the de voltage that is obtained after the pulsed voltage from windings ng and nu of the 
pulse transformer Tr2 is rectified by means of diodes D and Dy, and then smoothed out by capacitors Ci and Cz. The 
primary windings of output transformers Tr6 and Tr7 are connected to the collector circuits of triodes Ty and Ty. 
When triodes Ty, and Ty are unblocked, current pulses, which induce a voltage of complex form, us, at the outputs of 


Fig. 3. 
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transformers Tr6 and Tr7, pass through the windings nog and Ngg (Fig. 3e); after rectification, voltage uz is fed to the 
electronic keys for parameter switching (uy in Fig. 3f). 


The triode keys are based on an amplifier circuit with a common emitter, When negative voltage uy is fed to 
the bases of triodes Ty and ITI (Fig. 2, Ke-1), both triodes are unblocked, and the resistance for the large signal be- 
tween the triode collectors drops to 5-10 ohms. The other triode keys for the two current directions consist of triodes 
Ty and (Ke-2), Tg, and T (Ke~3), Tyg and (Ke- 5), and T and T (Ke-6). 


In most automatic control systems, the power consumption of the compensating network is low, This made it 
possible to construct a circuit based on low-power triodes P1-G and P2-B (Ty and Tyg), For protecting the input tri- 
odes from overloads, diode limiters (based on diodes Di, Dp, Dg, Dio, Dis, and Di and resistors Rg, Rg, Rig, RIS. Rag, 
and Rær) are connected to the inputs of the selective amplifiers; these diode limiters are supplied with pulsed blocking 
voltage from windings ng, Ng, and np. In this, the peak voltages that are fed to the collectors of all triodes do not ex- 
ceed 12 v. 


Fig. 4. General view of the control device. a) Magnetic ampli- 
fier; b) logical element unit. 


In our servosystem, the signals xi, xz. and xj, which are fed to the logical element unit, have polarities for 
which voltage uy appears in jumps if the conditions x} xj > 0 and x, x3 >0 are satisfied. 


Voltage x, is fed to the first input of the logical element unit and the input of the passive differentiating net- 
work, which consists of resistors B and Res and capacitor Cs. Voltage xz is fed to the second input of the unit 
(uinz = xz), and the signal of differentiation x, with variable parameters is fed to the magnetic amplifier. When volt- | 
age uy is fed to the input of Ke-1 or Ke · 2, the corresponding key opens, and the amplification factor K, increases in | 
Correspondence with Eqs. (10) or (11), The value of Ka, is determined by the magnitude of the variable resistance 
Rg. The voltage that is proportional to the velocity of the controlled variable %;, which is obtained from a tachom- 
eter generator, is fed to the passive differentiating network (Rx. Re and C.). 


The xj value, obtained by differentiation for constant parameters, is fed to the third input of the logical ele- 
ment unit (Ujns = x4). | 


Voltage x, is tapped off the unbalanced bridge- which is formed by fixed resistors Ry, Nur. Rs. and Ryg and the 
variable resistors Ry...» and Re-; of the keys Ke-2 and Ke-3~-and is then fed to the magnetic amplifier input. The 
arm resistance is chosen in such a manner that Ry + Rye-¢ Nes if there is no voltage at the key Ke-6 input. The 
magnitude and polarity of the magnetic amplifier input voltage are such that Ks K in correspondence with (17). 
If condition (15) is satisfied, the control voltage ug from transformer Tr6 is fed to the input of key Ke -G: however, 
since key Ke · I is blocked by the control voltage from transformer Tr7 for one current direction, the input of key 
Ke-6 is shunted, and the latter remains open (Ks K 
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If the sign relationship (16) is realized, key ke · I is deenergized, the shunting of Ke-6 is discontinued, voltage 
us from Tr6 closes key Ke · 6. and the ratio of the bridge arm resistances will be such that Ryg RC E- < Ri. As a 
result of this, the sign and the magnitude of the differentiating network output signal will change in accordance with 
(16). 


The error signal xi is fed to the input of the passive integrating network, which connects resistor Ras. capacitor 
Cs, and the unbalanced bridge, which consists of resistors Rai. Reg, Rag, and Rs. If condition (13) is satisfied, voltage 
uy Causes forced error integration as it arrives at the input of key Ke-2, since Ri + Rye-, < Rys RKe-s when key 
Ke-2 closes. At the same time, u arrives at the input of key Ke-4, which shunts the input of key Ke-3, and relation- 
ship (13) is realized independently of the sign of xz xj. 


For xi x C 0, key Ke-2 opens, and the integration is 
performed with the changed gain Kz Kg and time constant 
Tz (since the change in the time constant T2 is insignificant 
— in comparison with the gain Kg variation, the change in T, 
can be neglected in calculating the transient process), If 
relationship (14) holds, key Ke-4 is open, voltage uy from 
Tr6 closes key Ke-3, and the inequality for the bridge arms 
is reversed: Ryy + Rye-2 > Rus + RKe-3. In this case, integra- 
tion of the error with the opposite sign takes place accord- 
ing to (14), 


Figure 4 shows the general view of the control device- 
the summing magnetic amplifier and the logical element 
unit, This system was in operation under conditions where 
it tracked jump-like control inputs, which amounted to 6° 
on the final axle of the potentiometric bridge, 


Figure 5 shows the oscillograms of the servosystem 
transient processes and the diagram representing variations 
in the resistance of the keys connected to the circuit shown 
in Fig. 2 (Ke-3 and Ke-4), For condition (14), the resist- 
ance of key Ke-3 dropped stepwise almost to zero. In Fig. 5, 
curve 1 represents the transient process in a system with con- 
stant parameters (Ki 7, TI O. 15 sec, Kp=3, Iz = O. 6 sec, 
Ks= 0, Ks= 1; Tg= T;=0, K. 1. 17 * 1.9 sec, 16 * 0.1 sec); 
curve 2 represents the transient process in a system where 
che parameters change stepwise (K 10 3, KAy=2, Ty= 0. 15 sec, 
Keo=1, Kaz2=3, KA2g=9; the other parameters are the same); 
curve 3 indicates the variation in the resistance of keys Ke-3 
when the sign of the integrating element input quantity isre- 
versed, 


From an analysis of the oscillograms, it is obvious that 
the introduction of logical control according to the laws 
given by (9)-(14) makes it possible to obtain a system with- 
out considerable overcontrol, Although the introduction ofa 
passive integrating element does not make this system an 
¢ ‘ * astat ic system of the second order with respect to the control 

Fig. 6. Transient process curves, 1) Theoretical and input, the rather short integration time in the transient proc- 
2) experimental curves for a system vith constant pa- ess allows us to determine the basic dynamic properties of 
rameters; 3) theoretical and 4) experimental curves such systems, Due to the variability of the system param- 
for a system whose parameters change stepwise. eters, it was possible to reduce the transient process time 
sixfold. 


Figure 6 shows the oscillograms of the transient processes in a system with a second-order differentiating ele- 
ment (4 in Fig, 1) and the theoretical curves for the transient processes, which were plotted by using the Bashkirov 
graphoanalytical method. The parameter values for this system are the same as those obtained before. The param- 
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eters of the second-order differentiating element for curve 2 are: Kg =6, T}= 0.01 sec’, T}=0,007see; for curve 4, the 


parameters varied in the following manner: K 1. Kag=6. In this, the transient process time was reduced by a fac- 
tor of 4.5, and overcontrol was reduced sixfold. 


CONCLUSIONS 


1, We have presented one of the methods for the contactless variation of the compensating element param- 
eters in a servosystem by using a semiconductor logical unit with the aim of improving the transient process quality. 


2. The logical laws governing the variation of the system parameters make it possible to reduce four- to six- 
fold the transient process time and overcontrol. 
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ON THE THEORY OF SINGLE-TYPE INTERCOUPLED AUTOMATIC 
CONTROL SYSTEMS WITH SYMMETRIC CROSS-COUPLING NETWORKS 
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The principles of the synthesis of equivalent single-loop systems for single-type multichannel auto- 
matic control systems with symmetric cross-coupling networks are considered, The concept of syn- 
chronizing and averaging cross-coupling networks is introduced and a method for determining te- 
duced disturbances is given, 


Introduction 


Intercoupled automatic control systems with symmetric cross-coupling networks are used fairly often in prac- 
tice. Such systems consist of several identical automatic control channels which are intercoupled by means of equal 
cross-coupling networks. As an example, we can cite a system consisting of n identical electric power generators 
which are equipped with identical regulators and which operate on a common line; a two-channel servosystem for 
controlling the position of an axis in space; a system for controlling several duplicate · milling machines; computers, 
etc 
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The feasibility of developing a special theory for systems of this type, which would considerably simplify the 
stability analysis, was indicated in a paper by A. A. Krasovskii II]. This paper presented the concept of symmetric 
and antisymmetric cross-coupling networks in such systems and provided a method for investigating two-channel 
single-type systems with antisymmetric coupling networks by means of complex transfer functions. The present paper 
is concerned with the further development of methods for investigating such automatic control systems and with the 
extension of these methods to a special class of single-type multichannel systems which are characterized by the 
presence of the so-called synchronizing and averaging cross-coupling networks. 


The block diagram of a single-type multichannel automatic control system is shown in Fig. 1. In single - type 
automatic control systems, the input and output variables of identical elements are identical physical quantities, For 
the purpose of simplifying the analysis of such systems, it is convenient to introduce the concept of two types of mo- 
tion for the input and output variables of identical elements, The motion of the first type is the averaged motion, 
which correspond. to the average value of the input and output variables of identical elements, We shall assume that 
the number of channels in a single-type automatic control system is equal to n. We shall denote by i the number of 
the identical elements under consideration (in Fig. 1, the number of elements in a single channel is equal to m). The 
average value of the variable for the ith identical elements is equal to 


1 
a ru, 0 


where xi, is the average value of the variable of the ith element group, and xix is the variable value of the ith ele- 
ment in the kth channel. The average values of variables for all other elements are determined in a similar manner, 
Motions of the second type are relative motions which correspond to the difference between the input and output vari- 
ables of identical elements, These differences can be conveniently taken with respect to any single channel, for in- 
stance, the first channel, Thus, for the kth channel (Fig. 1), we can write 


= Ti · (2) 


It is obvious that for n channels there must be (n~1) such independent relative motions, Actually, a physical 
system whose parameters are strictly identical cannot be realized. Therefore, the question arises whether the ideali- 
zation which forms the basis of the consideration of single-type intercoupled systems is admissible: Is there a possi- 
bility that the small-value parameters which interfere with the identity of channels will cause substantial changes 
in the transient processes in the system ? 


t m 


Fig. 1. 


The small- value parameters which disturb the identity of channels in the system should be divided into two 
types. 


Small- value parameters of the first type (for instance, parasitic time constants) raise the order of the single- 
channel differential equation from N to N+v, l. e., they change the system's structure. In this case, the general dif 


; FD. 


sle- 
1 dif- 


ferential equation of the system will have additional roots, which are located at a considerable distance from the Nn 
basic roots, and the Nn basic roots are subject to small changes if the parameters of this type are small. When these 
mall - value parameters tend to zero, the basic Nn roots are shifted toward values which correspond to the identity of 
channels, and the additional roots tend to infinity. M. V. Meerov devoted a number of papers [2 and 3] to the investi- 
gation of the effect of small-value parameters of this type. M. V. Meerov found that the linear systems belonging to 
the large class which also comprises the multichannel systems considered in the present article are approximate in the 
sense defined by A. A. Andronov. 


Small- value parameters of the second type are such that they do not raise the order of the differential equation 
of a single channel, i. e., they do not affect the system's structure. Such small-value parameters comprise, for in- 
stance, small differences between gains or time constants in identical channels, The presence of small-value param- 
eters of this type changes the arrangement of the roots of the general differential equation; however, not all of the 
N(n-1) roots are multiple roots, It is known that the roots of a linear differential equation are continuous functions of 
its parameters [4}. Therefore, if the initial single-type automatic control system has a sufficiently large stability 
margin, small- value parameters of the second type cannot substantially change the character of the system's motion. 
Thus, the multichannel single-type linear systems considered in the present article are approximate in A. A. Andronov's 
sense, and a small departure from the identity of channels will cause a small deviation of the transient processes in 
comparison with the processes corresponding to the perfect identity of channels, 


Transfer Functions of the Equivalent Elements 


In order to investigate the averaged and the relative motions, it is necessary to compose equivalent block dia- 
grams for these motions. These block diagrams represent simple single-loop control systems, and they contain some 
equivalent elements whose transfer functions can be readily determined, The form of equivalent element transfer 
functions is determined by the character of the symmetric cross-coupling networks that are present in the system. 


Symmetric cross-coupling networks connecting identical elements are networks which transmit the signals from 
the outputs of identical elements to the inputs of the subsequent or preceding elements in other channels in the same 
proportion and with the same sign, Figure 1 shows examples of symmetric cross-coupling networks. There art forward 
and feedback symmetric cross-coupling networks. A forward symmetric cross-coupling network is a network through 
which the signal from the output of an element in one channel is transmitted to the outputs (or inputs) of the subse- 
quent identical elements in other channels. In Fig. 1, the first group of identical elements is provided with forward 
symmetric cross-coupling networks with the transfer functions LI (for the sake of clarity, Fig. 1 shows only three chan- 
nels), Symmetric cross-coupling feedback networks are those through which the signal from the output of an element 
in one channel is transmitted to the inputs of the preceding identical elements in other channels, In Fig. 1, the ith 
group of identical elements is provided with symmetric cross-coupling feedback networks with the transfer functions 
Lj. Let us first consider the case where the identical elements are provided with forward cross-coupling networks with 
the transfer functions L. From the definition of the forward cross-coupling network for the ith group of identical ele- 
ments, we can write the following n equations: | 


ix out = Wiexin + L > i in (k=1, 2. . ), 
1 


+k 


(3) 


where Xi} out is the output variable of the ith element in the kth channel, xiꝶ in is the input variable of the ith ele · 
ment in the kth channel, Wi is the transfer function of the ith element, and L is the transfer function of the forward 
cross-coupling network. By summing these n equations, we obtain 


> Zikout = >) Tin + 62 in — > uin). 
k=} k=] j=1 k=l 
Considering that 


n 
> Zixout = NA out, > Tixin = Mint, 
k=l k=1 


we finally obtain 
+ = [Wi + (n—1) in +. 
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Thus, in the presence of a forward symmetrical cross-coupling network for averaged motion, the transfer func- 
tion of the equivalent element in a single-type multichannel system is given by 


— 0 


For determining the relative transfer function of the equivalent element, we shall subtract the output variable 
of the row j (jn, j=#k) from the output variable of row in the system of equations GM. 


æu out out = M. (ruin in) + L (#45 in 5 


Hence, 


W,=L 


From the definition of the cross-coupling feedback network for the ith group of identical elements, we have 


By adding the left-hand and right - hand sides of the equations for all elements, we obtain 
n n n n 
Deu out = in + ru out) 
k=1 k=1 j=1 k=1 
or, by taking into account (1), 


WwW; 
Tou += 


Thus, the transfer function of the equivalent element for averaged motion in the presence of cross-coupling 
feedback networks is given by 


In order to find the transfer function of the equivalent element for relative motion, we shall subtract the ith 
row from the kth row in the system of equations (6): 


ru out — = Wi lzixin — zn L — tixout)]- 
Hence, the transfer function for relative motion is given by 


(8) 


Superimposed Cross-Coupling Networks 


The above expressions for the equivalent transfer functions for relative and averaged motion pertain to the case 
where identical elements are provided only with forward or only with feedback symmetric cross-coupling networks. 
In practice, cases are encountered where identical elements are provided with superimposed forward and feedback 
cross-coupling networks, 


Let us consider the typical cases of such superpositions of cross-coupling networks. Figure 2 shows the block 
diagram of a portion of a single-type multichannel system where the ith identical elements are provided with for- 
ward as well as feedback symmetric cross-coupling networks, In the block diagram of Fig. 2, the summation of sig- 
nals from the feedback cross-coupling networks is performed up to the point from which the signals to be transmitted 
to the inputs of the symmetric forward networks are tapped off, and the summation of signals from the forward cross 


— 


(4) 


(6) 


coupling networks is performed beyond the point from which the signals for the inputs of the symmetric feedback net- 
works are tapped off, For determining the transfer functions of the equivalent element for relative and averaged mo- 
tion, we shall write the equations relating the inputs and outputs of identical elements. For the ith identical ele- 


ments, we can write 
n n 
out WAu + Ly or out WA + LI — 
j=1 


(9) 
21. 2,..., a), 


where Aj, is the input variable of the ith identical element in the kth channel, and Li is the transfer function of the 
forward cross-coupling network. 
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For Aik, we can write 


= in + L. nj 
j=1 


‘ where I is the transfer function of the feedback cross-coupling network. 
Hence, 


Fig. 2. 


After adding the last expressions with respect to k (K 1. 2. 


n). 
we obtain 


(4 + W.L. IL) A = 2 ia 
1 


j=1 


(11) 


In order to determine the equivalent transfer function for relative motion, we shall subtract the row 4 from row 
k in the system of equations (9): 


ru out — Zijour = (Wi — Ly) Gu — 


By substituting the pertinent expressions (10) for Aik and Ajj, we obtain 


For determining the equivalent transfer functions for averaged motion, we shall add all n equations (9): 


k=1 k=1 
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n 
By substituting the found expression (11) for >) Ai, ve obtain 
k=1 


vl 
(13) 


2 rn in 
k=1 


A similar result is obtained for the second variant of coverage of the ith identical elements by superimposed 
cross-coupling networks, when the summation of the cross-coupling feedback network signals is performed beyond the 
point from which the signals of the forward cross-coupling networks are tapped off, and the summation of forward 
cross-coupling network signals is performed ahead of the point from which the feedback cross-coupling network sig- 
nals are tapped off, Another variant of superimposed forward and feedback cross-coupling networks is the variant 
where the forward cross-coupling networks are enclosed inside the feedback networks (Fig. 3). 


For the block diagram shown in Pig. 3, we can write 


Ain 
Tix out = = wd Tiʒ out + ny Tiʒ in + 
k 


+ (n— 2) LL + (n—1) LiLaxij out- 


| 


After e sat reduction of similar terms, we obtain 


(1 + — = (. Li) au in + 


+ (Wil, + (n — 2) LiLal Zijour + Ly in- (14 


j=1 


By subtracting the jth row from kth row, we obtain 
(1 + LiLi(æu out — Zijour = (Wi Li) (u in — 


Aout 


By adding all n —_ (14) and by reducing the similar terms, we have 
[1 —(n—1) Wily ru ot 
1 


Li >) tie in 


Ww W,+(n—1) 


— 

LESH 

—.ñĩ,é 


(14) 


1 in) 


(15) 


(16) 


By using a similar procedure, we can readily find the equivalent transfer functions for the last variant of super- 
imposed cross-coupling networks, where the feedback cross-coupling networks are internal and the forward ones are 
external networks. In this case, the equivalent transfer function for relative motion is given by 


1. . III. 


(17) 


while, for averaged motion, we have 


W, TGA - 1 WL, 
WL (18) 


Wi = 


Synchronizing and Averaging Cross-Coupling Networks 


In practice, cross-coupling between identical channels in single-type control systems is often a consequence of 
the difference between the output variables of identical elements, Thus, for instance, in a system of synchronous gen- 
erators that operate in parallel, the cross-coupling of generators through the common load appears only in the case 
where the rotor of one of the generators starts to run out (to lead or to lag) with respect to the rotors of the other ma- 
chines, Such cross-coupling between identical elements of single-type multichannel automatic control systems will 
be termed synchronizing cross-coupling [5]. Let us find the expressions for the transfer functions of equivalent ele- 
ments for averaged and relative motion for the case where the ith identical elements are encompassed by forward syn- 
chronizing cross-coupling. Figure 4 shows a section of a multichannel single-type control system with forward syn- 
chronizing cross-coupling networks (the synchronizing cross-coupling with the unit transfer function is conventionally 


denoted by a circle with a minus sign). By using the definition of synchronizing cross-coupling, we can write the fol- 
lowing n equations for the block diagram shown in Fig. 4: 


ru out = Wizwin + L >) (ruin (k= 1, 2,..., 


where W, is the ith element transfer function, and L is the synchronizing cross-coupling transfer function, This system 
of equations can be rewritten in the following manner: 


Zixout = [Wi (n — 1) — L > ij in · 


out 


Fig. 4. Fig. 5. 


W, 
Fig, 6. 


This system of equations describes the ith elements of a single-type system with ordinary forward symmetric 
cross-coupling networks, but with a different transfer function of the ith element [instead of M. the element transfer 
function is Wj +(n-1)L}. By using the above equations (4) and (5), we shall find the transfer functions for the equiva- 
lent elements: for averaged motion, 


Wiz = (Wi (n—1)(—L) =. (20) 


2 
— 
yin i — oy 


for relative motion, 
Wi = (Wy + (n— 1) L] —(—L) = W, 4+ nL. (21) 


On the basis of this, we can deduce that synchronizing cross-coupling does not change the character of averaged 
motion, since the equivalent transfer function for averaged motion is equal to the transfer function of an element ina 
separate channel, We shall find the equivalent transfer functions for the averaged and relative motion elements for 
the case where the ith elements are encompassed by synchronizing feedback cross-coupling. Figure 5 shows a section 
of a single-type multichannel system with synchronizing feedback cross-coupling. For such a block diagram, we can 


Zikout = 1. +L out 
1 


+k 


These equations can be rewritten in a somewhat different form: 


Ww n 


+k 


The obtained system of equations describes the ith elements of a single-type control system with ordinary feed- 
back cross-coupling networks, but with a different transfer function for the ith element. By using the above equations 
(7) and (8), we shall find the transfer function of the equivalent elements: for averaged motion, 


for relative motion, 
Ww, Ww; 


The same result has been obtained: Synchronizing cross-coupling does not affect averaged motion. 


In practice, we can also encounter cross-coupling which is the result of action of the sum of all output and in- 
put variables of identical elements, As an example, we can again cite synchronous generators operating in parallel, 
when their frequency (or voltage) regulators are connected to common bus-bars, In this, the frequency (or voltage) 
compensation will be provided with respect to the average frequency (or voltage) value for all the generators operat- 
ing in parallel, This type of cross-coupling between identical elements will be called averaging cross-coupling. 
Figure 6 shows the section of a single-type system with averaging forward cross-coupling (the averaging cross-cou- 
pling with the unit transfer function is denoted by a circle with a plus sign). From the definition of averaging cross- 
coupling, we have. 


1 n 
ru out in t+ 


This system of equations can be rewritten thus: 


This system of equations describes the ith elements of a single-type control system with an ordinary forward 
cross-coupling network, but with a different transfer function of the ith element IM + (1/n)L instead of WI. Accord- 
ing to (5) and (4), the equivalent transfer functions of an element with forward averaging cross-coupling are given by: 


—— 


(23) 


(24) 


(25) 


for relative motion, 

L L 
Wy = Ws 
for averaged motion, 


Wy TG =. I. 


It is obvious that averaging cross-coupling does not affect relative motion, By using a similar procedure, we 
can readily find the equivalent transfer functions for the case where the identical elements are encompassed by aver- 
aging feedback cross-coupling. For this, the equivalent transfer function for relative motion will be 


W. 


(28) 


while, for averaged motion, we shall have 


— - 
(29) 


A similar result is also obtained in this case: Averaging cross-coupling does not affect relative motion, 


It can happen that, along with simple symmetric cross-coupling, synchronizing and averaging cross-coupling is 
also present in a single-type intercoupled automatic control system. Considering the simplicity of investigating inter- 
coupled single-type systems where only synchronizing and averaging cross-coupling is present, it is convenient first to 
reduce all simple cross-coupling networks to synchronizing and averaging networks. Let us first consider how a simple 
forward symmetric cross-coupling network should be decomposed into a synchronizing and an averaging network. For 


the ith group of identical elements with forward symmetric cross-coupling, the n equations (3) hold, Let us transform 
the last term in this system of equations: 


n n n 
L > in = in + os 


{i j=1 


n n 
— L 
=* L 2 in — 22 (ru in in)- 
j=1 


The system of equations (3) can now be rewritten thus: 


15 15 
ik out Miru in + (n—1) — (ain in). 


If we compare the last equation with the systems (19) and (25), we notice that the simple forward symmetric cross- 


coupling network is now decomposed into a synchronizing network with the transfer function —L/n and an averaging 
network with the transfer function (n=1)L. 


A symmetric feedback cross-coupling network can be decomposed in a similar manner, For this, we shall trans- 
form the last term in the system of equations (6): 


LD Zour =(n—1)L+> >) Lruout = 
=1 j=1 j=l 


+k 
15 ix 
j=1 
The system of equations (6) can now be written thus: 
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Thus, the feedback cross-coupling network is now decomposed into a feedback synchronizing network with the 
transfer function Ln and feedback averaging network with the transfer function (n-1)L. 


Equivalent Disturbances 


In determining the transfer functions for averaged and relative motion, nothing has been said about disturbances 
that may be present at the inputs of identical elements, By using the superposition principle, we can formulate a rule 
for determining reduced disturbances for equivalent elements for averaged and relative motion. 


We shall first assume that the disturbance qj, acts at the input of the kth element in the ith group of identical 
elements, This disturbance can be represented in the form of the sum 


n—i 


We shall now assume that ident ical and opposed disturbances 2 Qj, act at the inputs of the identical elements 


with the exception of the kth element, Then, we can consider that two disturbance systems act at the identical ele- 
ment inputs, The first disturbance system consists of n equal disturbance functions (/n) dig -· This disturbance system 
must be taken into account only in investigating averaged motion; it does not affect to any extent the relative mo- 
tions, Thus, the equivalent disturbance for averaged motion will be (17 n) qi · The second disturbance system consists 
of n- equal disturbances ~(1/n)qj,, which act at the inputs of all elements with the exception of the kth element, 


and of the disturbance — Gik, which acts at the kth element input. The sum of all these disturbances is equal to 


zero, and, therefore, the latter disturbance system does not exert any influence on averaged motion. For all relative 
motions with the exception of the kth motion, the equivalent disturbance for the second disturbance system is equal 
to zero, and, for the kth relative motion, the equivalent disturbance is equal to 


1 2 —1 


Following a similar procedure, it can be readily shown that, in the general case, when disturbances act at the 
inputs of all identical elements, the equivalent disturbance for averaged motion is equal to 


1 n 
Tint = = >) 
j=1 


while, for relative motions, we have 
Jij— = Vi (33) 


In investigating the quality of control processes for relative motion, the greatest of all the 4n(n-1) disturb- 
ances should be selected, since this would represent the most adverse case, 


As an illustration of the practical application of the concept of synchronizing and averaging cross-coupling, we 
shall consider here an example of the composition of equivalent block diagrams for the averaged and relative motions 
of a single-type control system consisting of n identical synchronous generators which operate in parallel with a com- 
mon load and which are equipped with speed regulators, excitation regulators, and devices for the uniform distribution 
of resistive and reactive power. The schematic diagram for the parallel operation of single-type synchronous genera- 
tors is shown in Fig. 7 (for clarity, only two generators are shown), The synchronous generators G are driven by piston 
motors, turbines, or hydraulic drives, which are provided with speed regulators SR. The synchronous generators are 
equipped with excitation regulators ER, Facilities for the automatic equalization of the resistive and the reactive 
power, based on the virtual statism method, are provided in the system [6]. Frequency compensators FC, which act 
on the speed regulators of the generator drives, serve for the automatic resistive power distribution. Voltage com- 
pensators VC, which act on the generator excitation regulators, serve for the automatic reactive power distribution. 
The frequency and the voltage compensators are controlled by the resistive load distributing device ResLDD and the 
reactive load distributing device ReaLDD. The frequency and voltage compensators are also controlled with respect 
to the average frequency and voltage values in the common bus-bars, 
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The following assumptions were made in composing the system equations; 


the 

1, All of the n units that operate in parallel are identical. 

2, The generators operate with a fixed resistive load in the line, which is located symmetrically with respect to 
all the generators, 

pre 3. The effect of the stator aperiodic current components and the corresponding rotor periodic current compo- 
nents is not taken into account, 
ical 4, It is assumed that the generators are not provided with antihunt windings. 

5. The synchronous generators have equal reactances in the longitudinal and transverse directions (the resist- 
ance of the stator windings is neglected), The last three assumptions are not basic, and the equations can be written 
without these assumptions, 

The equations for the kth unit in operator form for the zero initial conditions are given below. The equations 

nents are written in relative units, while the steady-state values of the variables are taken as the basic quantities, 
ele- The drive equation for the kth unit is 
ystem 1 n n 
j=1 j=1 
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The equation of the speed regulator for the kth unit is 
ive 
h act % = N, (8) ( + Cy): 
om- 
ton. The equation of the frequency compensator for the kth unit is 
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The equation of the electromagnetic précesses in the kth unit generator is 


n 
The equation of the excitation regulator for the kth unit is 


= Nu (e) + 


The equation of the voltage compensator for the kth unit is 


nn n n 

ugjt Le > (4, — + Lio (v, — ¥;)- 
j= 
i+k 


In all these equations, vy, is the relative rotating speed of the kth unit, Oj, is the relative output variable ofthe 
kth unit speed regulator, & is the relative output variable of the kth unit frequency compensator, uc, is the relative 
voltage at the terminals of the kth synchronous generator, edek is the relative output variable of the kth unit excita- 
tion regulator, E is the relative output variable of the kth unit voltage compensator, q; is the relative disturbance, 
A, Ly, Lg, Lg, Ly, Lg, Lg, Ly, LS. Lo. and Lie are coefficients (generally operator coefficients) which depend on the sys- 
tem parameters and which are determined by analyzing the equivalent circuit and the vector diagrams of the synchro- 
nous machines operating in parallel, R,,(s) is the speed regulator transfer function, and R,,(s) is the excitation regu- 
lator transfer function, 


Fig. 8. 


Figure 8 shows the block diagram of a system which is based on the above equations, For the sake of simplic- 
ity, only two units are shown in this block diagram. It follows from the above block diagram that a system consisting 
of n identical synchronous generators which operate in parallel represents a single-type multichannel automatic con- 
trol system with symmetric synchronizing and averaging cross-coupling networks, 

In investigating the stability of the system's relative motion, all averaging cross-couplings must be excluded 
from consideration, By using the detecting properties of the elements and the superposition principle, the block dia- 
gram for relative motion can be reduced to the form shown in Fig.9. By combining the drive element with the speed 
regulator element and the generator element with the excitation regulator element, we obtain a block diagram rep” 
resenting a series assembly of three groups of identical elements, 


15 
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The first group of identical elements with the transfer functions sT = R- 1+ Ry (s)]"* is encompassed by nega- 
tive synchronizing feedback couplings with the transfer functions Lg+LsR,(s). The second group of identical elements 
with the transfer functions A is encompassed by negative forward synchronizing couplings with the transfer functions 
igRy(s). The third group of identical elements with the transfer functions Is I 1+ Ry () l is encompassed by nega- 
tive synchronizing feedback couplings with the transfer functions L Ru (). All of the three groups of identical ele- 
ments are encompassed by negative synchronizing feedback couplings with the transfer functions LV. Liu (s). 


Fig. 10. 


TI f . By using the above equations (21) and (24), ve find 
5 Ra the transfer functions of the individual groups of elements 
and of the system as a whole, The equivalent transfer func- 
— | tion of the first group of identical elements is given by 


Ww, = = 
+ LR, ler. + R, 
= [sT, N, (8) (4 + nls) + 
The equivalent transfer function of the second group of identical elements is 


Wy = A+ LR, (8). 


The equivalent transfer function of the third group of identical elements is 


[1 Ry 
[1 + (s)] [1 + Ru = 1+ Ry % (1 + 


W = 


The equivalent transfer function of the entire closed-loop system for relative motion is given by 


W Wy Wn A+ (8) 


Sy = [sTy + R—-A+R, (s) (1 + + [1 Ru (1 
S. [Ly + (8)] [A ALR. (l. 


Thus, the conditions for stable relative motion of a system consisting of n single-type synchronous generators 
which operate in parallel can be reduced to the requirement that the real parts of the roots of the following charac~- 
teristic equation be negative: 


+ R—1+R, % (1 + Ry (9) (1 + + 
+ n[L, + LioR, (s)] [A + LN, = 0. 


The conditions for stable averaged motion can be determined in a similar manner, In this, the synchronizing 
ctoss-couplings must be excluded from consideration, The block diagram for averaged motion is shown in Fig. 10. 
The equivalent transfer functions for all of the three groups of identical elements can be found by using Eqs. (27) and 
(29), The equivalent transfer function for the first group of identical elements is 


WI. =([sT, + R—1+ RN, 


| 


The equivalent transfer function for the second group of identical elements is 
A — I- 


The equivalent transfer function for the third group of identical elements is 
= 1 IL. + Ry (8) (1 + 


Since the considered system of three groups of identical elements for averaged motion represents an open-loop 
system, the averaged motion stability conditions can be reduced to the requirement that the real parts of the roots of 
the follow ing characteristic equations be negative: 


R- R. =. +1 L. — N. (8) (1 + Le) =0. 


It can be shown that these conditions are adequate as stability conditions for the speed and voltage control sys- 
tems of a single unit under the assumption that the sensing elements of the speed and voltage compensators are con- 
nected to the bus-bars proper of the generator of a single unit. 


CONCLUSIONS 


The analysis of single-type multichannel automatic control systems can be considerably simplified by reducing 
them to equivalent single - loop systems. A method for composing such equivalent single - loop systems has been pre- 
sented in this article. The described method for determining equivalent disturbances makes it possible to investigate 
not only the stability, but also the performance of automatic control systems by using the above single-loop systems, 
The concept of synchronizing and averaging cross-coupling networks has been introduced, and a rule for decomposing 
simple cross-coupling networks into synchronizing and averaging networks has been presented, 


The investigation of single-type intercoupled automatic control systems with synchronizing and averaging cros- 
couplings has been greatly simplified due to the fact that relative motion is invariant with respect to signals that have 
passed through averaging cross-coupling networks and with respect to disturbances, which are equal for all channels, 
while averaged motion is invariant with respect to signals that have passed through synchronizing cross-coupling net- 
works and with respect to disturbances, whose sum is equal to zero for all channels, 
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THE TRANSFER FUNCTION OF AN AUTOMATIC CONTROL SYSTEM 
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An automatic control system with a modulator and a half-wave demodulator is reduced to an equiva- 
lent continuous-control linear system. The method of modulated harmonics is used for determining 
the transfer functions, 

The effect of converter elements consisting of a modulator and a half-wave demodulator on the 


dynamics of a servosystem with a crossover frequency of 25cps is investigated, The results obtained 
in the experimental verification of the theoretical data are given. 


1. Introduction 


Consider the automatic control system whose block diagram is given in Fig. 1. The modulator M,. the half- 
wave demodulator DM, and the linear element L form a closed group, Disturbance xg acts on the system. 


The harmonic signal g with the carrier frequency w», which is modulated by signal @ and is fed to the modu- 
lator input, is produced at the modulator output: 


6 (t) = K. t) sin apt. (1) 


The demodulator rectifies the signal $(t) with respect to the sign of the modulating signal t). It is assumed 
that the coupling between the modulator and the demodulator is inertialess, 


The demodulator output signal is expressed by the equation 


1 (t) Kit) A (ont). (2) 
where 


A (Wot) = sin dt for 2nn < ag <(2n + IA. 
(n =0, 1. 2, 3) . 
A (Wot) =0 for (2n +1) < aot En + 2) x. @) 


The demodulator output signal n (t) consists of a sequence of positive sinusoidal half-waves, which are modu- 
lated by the mismatching ot). The diagrams of ot), A (t) and n(t) are given in Fig. 2. 


We shall assume that the transfer function 


7 (s) 0 


for the element L is known, 


The system under consideration belongs to the class of unste ady- state linear systems with a periodically chang- 
ing parameter. Such systems are often encountered in automatic control technology. As an example, we shall cite a 


servosystem where selsyns play the role of the summing unit and the modulator, a diode phase detector acts as the 
demodulator, and an amplifier and a motor represent the linear element, 


The presence of a section with a modulator and a demodulator for which Eq, (2) holds in the automatic control 
loop entails certain difticulties in determining the transfer function of the closed-loop system, In analysis, it is cus- 
tomary to circumvent these difficulties by replacing the half-wave demodulator by an ideal one, In using this ap- 
proach, the modulator input quantity is related to the demodulator output quantity by the equation 


(t) = (i). (5) 
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where Kp are the amplification factors of the modulator— demodulator element. 


The present article is concerned with the determination of the closed-loop transfer function where the relation- 
ship between the modulator input quantity and the demodulator output quantity as expressed by Eq. (2) is taken into 
account and also with the determination of the conditions under which the replacement of the half-wave demodula- 
tor by an ideal one is admissible. 


We shall solve the stated problem by using the method of modulated harmonics [1], which we have used before 
for determining the transfer function of a step automatic control system. 


2. Expansion of the Input and Output Quantities of Ele- 
ments into Modulated Harmonic Series 


Tq 

9 As was mentioned above, the demodulator output signal ij (t) represents a 
sequence of positive half-waves of the A (t) sinusoid, which is modulated by 
mismatching et). The function A (tt) can be expanded into a Fourier series 
and represented as a harmonic sum. The demodulator output signal ij (t) can 
then be considered as a harmonic series where the harmonic amplitudes varyin 
time, 


The 7 (t) harmonics will pass through the linear element L and the sum- 
ming unit, The output signals x(t) of the linear element and @(t) of the sum- 
ming unit will also consist of harmonics whose amplitudes change in time: 

+00 
y(t) = > N. (t) 
k=—0o 


= 
wt ＋ (t) * x (t) 2 


(t) = (i) 
k=—oo 


rg, 1 Time diagrams of the mis- 
tt fi t f 
matching @, — fenc- The nu (t). xx (t). and e (t) functions define the variation in time of the 


signal. We shall assume that the functions nu (t), xh (t), and & (t) for all k from 
— to satisfy the conditions under which Fourier transformation is possible. 
We shall denote by u (O). xk (), G () and xqg(w) the Fourier images of the functions nx (t). xK(t), O(t), and 
xd (t). 


The ratio of the image of the controlled variable kth harmonic amplitude to the disturbance image, 


() 


Y.(@)x = 


will be called the complex transfer function of the closed-loop system with respect to the kth harmonic. 
The problem consists in determining the transfer functions Ie (w),. 


We shall write the equations relating the images of the amplitudes of harmonics , x, and @ to the image of 
Xq- 


The summing unit equations will relate x,(w) and 6 (w) to xꝗ (o). the modulator— demodulator equations will 
relate q (w) to n (w), and the linear element equations will relate Ny, (w) to x, (). 


a) Equations of the Summing Unit 
Let us substitute in the summing unit equation 


(t) + 6 (t) = za(t) 


A 
+0o 
(6) 
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(6) 


the values of x(t) and ot) from Eq. (6): 


2 lar (t) + % (t)] 24 (f). (9) 


Equation (9) will be identically satisfied for any t if the amplitudes of the harmonics are balanced: 


Zo (t) + (%) = (t), (t) +O (t)=0 (#0), (10) 


The balance equations for the harmonic amplitude originals will correspond to the balance equations for the 
images: 


Tq (o) + (w) = (w), (o) +O () = 0 (11) 
b) Modulator Demodulator Equations 


Let us expand the function A (ot). which is expressed by Eq. (3), in a Fourier series that is given by 


+co 
A (wot) = >) 


1 
0, 4 2, +4,...; 


c for n= 43, +5, +7,...; 
e 


(13) 


By substituting @(t) from (6) and A(w¢t) from (12) in Ed. (2), we obtain 
+00 +00 


After performing the multiplication operations and grouping the terms with equal exponents, we obtain, in cor- 
respondence with (6), the equations relating the amplitudes of harmonics n and 6: 


eo 
K. e. b. (k = — to. (15) 


By applying Fourier transformation to (15), we obtain a system of equations which relate the images of the am- 
plitudes of harmonics 9 and @: 


(@) = K, > cn ( (k= — oo), (16) 


c) Linear Element Equation 


Let us denote by nu, nt) and Xk, wa (t) the kth modulated harmonics at the linear element input and output: 


n. „(t) = 1, (t) e, „(t) = (t) e. 1 
The images of the modulated harmonics will be denoted by n. w.(w) and x, .,,(w), respectively. 


Let us apply Fourier transformation to Eqs, (17),denoting the variable in the region of images by My. By taking 
into account the the displacement theorem, we obtain 


Nk (2x) (2. k@»), (Q.) Zp (Qk — k@o). (18) 
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The images of the linear element input and output quantities are related by the linear element complex trans 
fer function V. (an): 


(2,) VI. (Q,) Nk. „). 


(19) 


From (18) and (19), we obtain 


x, (Q, — = Ii (2,) n, (Q. — ke). 
Wie substitute for the independent variable of the image region by putting 


— = . 


Then, 


x, (o) = II (@ + (o) ( = — 40 05 


The system of equations (22) establishes a relationship between the images of the harmonic amplitudes at the 
linear element output and input. 


‘ 
3. Transfer Functions for the Amplitudes of the Controlled Variable Hatmonics 


Equations (11), (16), and (22) relate the images of the harmonic amplitudes of signals t), n (t), and x(t) tothe 
image of disturbance xq(t). 


By eliminating the images of the amplitudes of harmonics n (w) and d (), we obtain a system of equations 
which relates x; (w) to xq(w): 


+00 
(h + (@ + Cnn = (@ + hog) 24 (o) 


n=—oo 


(k = — o to + oo). (23) 


The system of equations (23) can be solved for any xx, and, in correspondence with (7), the closed-loop system 
transfer functions for the harmonic amplitudes are found. 


In the general case, the system (23) consists of an infinite number of equations. In most automatic control sys 
tems that are encountered in practice, the linear element L has the properties of a low-frequency filter, and it causes 
an intensive attenuation of the input signal components, which have relatively high frequencies, In such cases, xy) 
and M. (w+ kw»), which correspond to higher harmonics, can be omitted from the system of equations (23), We shall 
consider two particular cases, 


a) Only the Zero-Harmonic Amplitude Is Retained in the System of Equations (23) 
In a certain frequency range, let 


I. (@+ ke) (k*0), 


The system of equations (23) then assumes the following form: 


25 (o) + (0) 4, (0) Kier. (o) 24(0), * 


Klee (o) = I (o) (= (oh, 


where 1 (o) is the open-loop linear system transfer function, and Kp is gain of M and DM with respect to the direct 
component. 


|_| 
(21) 
— 
(24) 
(26) 
Let us denote 
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(23) 
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From Eqs. (7) and (25), by taking into account (27), we find the closed-loop system transfer function for the 
zero-harmonic amplitude: 


10 (w) Y (@) 
(@) — 44 (0 — 1+Y(@) . 


(28) 


Expression (28) corresponds to the transfer function of a linear closed-loop continuous-control system which con- 
sists of element L with the transfer function II (C) and an inertialess element with the transfer function Kp. which re- 


places the modulator and the demodulator. Equation (28) corresponds to the case of an ideal demodulator IK C). 


The stability of a closed-loop system with the transfer function (28) can be estimated with respect to the trans- 
fer locus Y (H) of the open-loop system by using the Nyquist— Mikhailov criterion. 


For the determination of stability, it is especially important to know positively the transfer locus portion which 
lies close to the critical point whose coordinates are —1, jo. For the crossover frequency w=Weo, the transfer func- 
tion modulus is equal to unity, The transfer locus in the crossover vicinity will be determined in a sufficiently reli- 
able manner if the following conditions are satisfied: 


Y (@ggt (k 0). (29) 


For a positive Weo value and for k=—1, the most adverse condition will be 


Y (@eg— @o) 0. (30) 


In calculating the closed-loop system transfer function for the zero-harmonic amplitude, the higher harmonics 


may be neglected if, for the Weg- we frequency, the modulus of the open-loop system transfer function is small in 
comparison with unity. 


b) Only the Amplitudes of the Zero and the First Harmonics Are Retained in the System of Equations (23) 


In a certain frequency range, let 
Y(@+ko@)=0 


The system of equations (23) then assumes the following form: 


(31) 


[co + K (e) + (0) + en (0) = 
(@) + [co + | (@) + cri (@) = (). (32) 


+ (0) + [60+ (0) = 


We shall solve the system of equations (32) with respect to x9(w), x-,(w), and xi (). and, by using relations 
(13) and (27) and taking into account (7), we shall find the closed-loop system transfer functions for the zero and the 
first harmonic amplitudes: 


¥ () 


(@) Yes, 
YX (oa = T „Je (% = Te) 
(35) 


2 
Ve ( + @o) + Ve (@—@o) J (@+ ( 
(w) = ¥ (o) [1—(4) (36) 
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1 
1— V (@ + 0 
(o = le (% 1 2 
1—% (@ + @) (% — 


1— 


Yo(w 


where 1 (o) is the open-loop linear system transfer function which is expressed by Eq. (27), and Yo (w + wg) and 
— are the closed-loop linear system transfer functions for the w+ Ww» and w —w» frequencies, respectively, 


Yes (@) = le (o 


Y (o Y (@ — o) 

Yo(@ + = Je (o = TTT (38) 

The closed-loop system transfer function Yo (w), for the controlled variable zero harmonic, which is expressed 

by Eq. (34), corresponds to a block diagram where the element with the transfer function Ye (w) is provided with nega- 
tive feedback (Fig. 3a). According to Eqs, (35), the closed-loop system transfer functions for the harmonics with k==1 
and k= 1 correspond to the block diagrams shown in Figs. 3b and3c, Thus, a closed loop consisting of a modulator, a 
half-wave demodulator, a linear element, and a summing unit is equivalent to a linear closed loop without the modu- 
lator and the demodulator. The transfer function Ye (w) of the equivalent open-loop linear system is expressed by 
Eq. (33), and the transfer functions of the equivalent elements whose output quantities correspond to the amplitudes of 
harmonics with k= 1 and k=1 are expressed by Eqs. (37). 


Yo) Selo) 
b c 


Fig. 3. Block diagrams for the zero (a) and the first ( and c) har- 
monics of the controlled variable. 


It is obvious from Eq. (36) that Ye () differs from the open-loop linear system transfer function Y(w), which 
corresponds to the case with an ideal demodulator. 


Let us further assume that, for positive w values in the frequency range wS Weg (where co is the crossover fre- 
quency for which Iv (col = 1), the following inequality is satisfied: 


| ¥ @) |< 1. (39) 
Then, according to (38) and (39), 
Te (Wegt | <1. * 


and the terms containing Yc (w+ ) can be omitted in (36) and (37). In this case, the transfer functions of the equiv- 
alent linear elements will be expressed by the following equations: 


Ye (o) = (o) — (D' Yo (o (41) 


Te () = Tete — a), Yer (o) = 0. (42) 


The stability of a closed-loop system with a modulator and a half-wave demodulator can be determined with 
respect to the transfer locus Ye (o) for the open-loop linear system by using the Nyquist— Mikhailov criterion. 


It follows from (41) that the larger the modulus of the expression 


(Tre (@co- oo). 


— 


the greater the effect of pulsations at the demodulator output on the transfer locus in the vicinity of woo and, conse - 
quently, on stability. 


4, Investigation of Servosystems with a Half-Wave Demodulator 


It was shown above that, under certain conditions, the transient processes can be considerably affected if the 
automatic control system contains a converter element consisting of a modulator and a half-wave demodulator which 
is not an ordinary linear continuous-action element. 


The obtained results were checked experimentally by means of a servosystem which consisted of a selsyn meas- 
uring element, a half-wave demodulator, an electronic amplifier, and a de motor. 


The system performance was investigated for the 400-cps and 50-cps supply voltage frequencies of the selsyns. 
The system parameters were chosen in such a manner that the moduli of the transfer functions Yo (weg + wg) and 
Yo(Weo ~ had a value much smaller than unity for the frequency fy=400cps and a value of the order of unity for 
fy=50cps, The time constants and the gains of the system elements were equal for both frequencies fo. 


Servosystem Wiring Diagram (Fig. 4) 


The servosystem measuring element is based on selsyns, which operate in conjunction with transformers, The 
output voltage of the selsyns is fed to the half-wave demodulator, which consists of transformers TI and T2, a twin 
diode, LI. and load resistors Rg and Rg. The resistor Ri in the T1 primary winding circuit serves for regulating the de- 
modulator gain, An ac reference voltage is fed to the T2 primary winding, The demodulator output voltage U, is 
added to the flexible feedback voltage Us (which is provided by resistors Rig) and then fed to the cathode follower L2. 
The cathode follower output voltage is fed to a three-stage push · pull electronic amplifier, whose preamplification 
stages are based on tubes L3 and LA and whose end stage is based on tubes L5-L8, Filters Cg, Ry, and Rg, which serve 
for smoothing out the pulsations of the signal arriving at the amplifier input, are connected to the grid circuits of the 
L3 arms, Potentiometer Ry serves for balancing the amplifier. 


The amplifier output voltage Us is fed to the final separate-excitation motor Mo. The motor armature is cou- 
pled through reducing ge ar Re with the rotor of the receiver selsyn RS and the armature of tachometer generator TG. 
Voltage Uy from the TG armature is fed to the differentiating feedback circuit CRs, which is used for stabilization, 


The servosystem block diagram is shown in Fig. 5. The servosystem elements have the following units: Ele- 
ment 1 includes the selsyns and the demodulator; element 2 consists of the electronic amplifier; element 3 combines 
the motor and the reducing gear (the output signal of element 3 is the reducing gear angular velocity w); element 5 
consists of the tachometer generator; element 6 consists of the differentiating circuit CRE. 


Element 4 acts as an integrator, and it converts the reducing gear shaft angular velocity into the rotation an- 
gle ?. 


The transfer functions of the elements are given by 


Y, (@) = K, 
(for the direct component), 


(10 
7. (0) = Yao) 
(43) 
r. (ig, K., V. (eh 


The experimentally found gain values are equal to: Ki = 0,36 v/deg, Kg = 900, Kg = 118 deg/vsec, Kg = 
= 1,8-10°* vsec/deg, 


The feedback element time censtant is T = RISC = 17.5+ 107? sec, 
The amplifier filter time constant is T, = NC 17.5 · 10˙ sec; T, = “aS = 0,106, 


The motor time constant Tm was found with respect to the experimentally determined frequency response of 


the circuit portion consisting of the terminal amplifier stage, the motor, and the tachometer generator; it was Tm = 
= 0,106 sec, 
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Transfer Functions and the Logarithmic Amplitude and Phase 
Responses of the Servosystem 


Elements 2 and 3 (Fig. 5), which are provided with feedback 
through elements 5 and 6, can be represented in the block diagram 
by a single equivalent element with the transfer function 


Y2 (@) Ys (o) 
Y,(@) = 1+ Ys 00 Ys ( Ve () (44) 


By substituting expression (43) in (44) and considering that the filter time constants are TI Tm and T = T. ve 
obtain 


uber +1) 


It is obvious that the equivalent element is stable. Let us first consider the servosystem as a continuous · contiol 
system while neglecting the effect of pulsations of the converter element output signal. 
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The open-loop system transfer function will be equal to 


¥ (@) = (0)¥s (0) K 


Here, K = Kikzk, = 38 · 10° sec” is the open-loop system gain. 


Fig. 6. 


The open-loop system logarithmic amplitude and phase responses are given in Fig. 6 (curves Y). It follows from 
these curves that, for the crossover frequency fgg= 25 cps, a phase margin Ag = 46° is available, which secures ade - 
quate damping of the closed-loop system free oscillations [2]. 


If the selsyns are fed from a line with the frequency fy = 400cps, then fog —fy =—375cps, and fog+ % 425 cps, 
For these frequencies, the open-loop transfer function is considerably less than unity, and the transfer function Ie (w) 
of the equivalent open-loop system degenerates into Y(w), For the line frequency f,=400cps, the servosystem will 
behave as a continuous-control system with the frequency responses marked by Y in Fig. 6. 


Experimental investigations showed that the system is stable for fg=400cps, An oscillogram of the tachometer 
generator voltage when the servosystem processes a jump-like mismatching of @°=5° is shown in Fig. 7. 


If the selsyns are fed from a line with the frequency fg=50cps, foo — fp = 25 cps, and fog+ % 75 pe. As fol- 
lows from the Y curves in Fig. 6, the modulus of Y(w) in the vicinity of foo fy is close to unity, and, consequently, 


the transfer function Ye (w) of the equivalent open-loop system will differ from the transfer function Y (w) of the lin- 
eat open-loop system. 
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Figure 6 shows the logarithmic plots of the amplitude and phase responses that correspond to the transfer func- 
tion Ye (w) of the equivalent open-loop system, which was calculated by using Eqs. (38) and (36) (curves Ye). The 
same figure shows the response curves corresponding to the simplified expression for the transfer function Ig (), which 
was calculated by using Eqs. (38) and (41), It follows from these graphs that the equivalent system will be unstable, 
In the vicinity of the crossover frequency, the le ( and Yg(w) curves are rather close to each other. 


rig. 7. 
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Fig. 8. 


Experimental investigations confirmed the instability of the system. Figure 8 shows an oscillogram of the tach- 
ometer generator voltage. It is obvious from the oscillogram that, for fy=50cps, undamped oscillations with a 25-cps 
frequency are present in the system. The determination of the frequency and amplitude of undamped oscillations isa 
nonlinear problem which does not fall within the scope of the present paper. 


5. Conclusion 


The input and output quantities of each element in an automatic control system with a modulator and a half- 
wave demodulator can be represented in the form of a harmonic series where the harmonic amplitudes vary in time. 
The amplitude of each of the harmonics corresponds to an equivalent closed-loop system consisting of continuous lin- 
ear elements. The complex transfer function of the equivalent open-loop system is expressed by the transfer function 
Y(w) of the open-loop linear system and the closed-loop linear transfer functions I (W + kw») with the har- 
monic side frequencies. 


It the crossover frequency of a linear open-loop system is much lower than the carrier frequency, then 
I Weg + kw») « 1 for k= 0, and the equivalent system transfer function Ye (o) degenerates into the linear system 
transfer function Y(w), In this case, the demodulator pulsations hardly affect the system behavior, and the demodu- 
lator can be considered as an ideal demodulator. 


If the modulus | Yo (c k is commensurable with unity, if only for k =—1, the equivalent system trans- 
fer function will be different from Y(w), 


The theoretically predicted influence of a converter element consisting of a modulator and a half-wave demod- 
ulator on the transfer function of the system was qualitatively confirmed by experimental investigations, In the case 
under consideration, the system, which was satisfactorily stabilized for the modulator carrier frequency o- 400 cps, 
became unstable for fy=50cps, 
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A GENERALIZED GRAPHICAL METHOD FOR DESIGNING 
CONTACT CIRCUITS 
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A generalized graphical method for designing contact (p, k)-poles by taking into account bypass paths 
is described. The procedure to be used in formulating the conditions for the ambiguous solution cases 
and a method for selecting the optimum solution are given, 


1. Introduction 


Along with the algebraic methods for designing contact circuits, “direct” methods are also being presently de- 


veloped; from among these methods, we can mention general graphical methods [1-4] as well as methods for the par- 
ticular case of symmetrical circuits [1 and 5]. 


For designing a (1,k)-pole circuit by using the graphical method, all networks are given as number sets 
fx = Nas - „ Nr, (ur, 


(1) 
of those states for which a given network must be (“obligatory” numbers n) or may be ("conditional” numbers n) 
closed from the output x to the common input. By the state number, we understand the sum of "weights" di = 211 
of those relays which operate in a certain given state (i=1, 2, 3. . n is the relay index number). The numbers 
belonging to the set also depend on the circuit relay numeration order (assignment of weights), For n relays, there 
can be nf such orders, The chosen order will be called the “base” B and noted in the form of the set index, while the 
relays in the base will be arranged in the order of their diminishing weights, Thus, B = XXX, signifies that weight 1 


has been assigned to relay XI, weight 2 to Xa, and weight 4 to Xs, In order to simplify the notation, we shall denote 
the base X; . . . XI simply by i. 


Expression (1) indicates that the given conditions can be satisfied by 2° different contact circuits, which differ 
from each other by the combination of the numbers of conditional states in which the circuit will be closed, Thus, 
the structural admittance ¶ of the circuit to be designed must lie within the interval 


A set which includes the weight of 2! numbers (for the I base) will be called a complete set, while, if the con- 
ditional numbers are not present in this set, such a set will correspond to a permanently closed circuit, i. e., {o, 1, 
eses 2—i}is 1, A set without compulsory numbers will be called an empty set, In the case where even the con- 
ditional numbers are absent, the set will correspond to a permanently open circuit, i. e., {-}= 0, 


The sets where all compulsory numbers of any set are contained in each of the remaining ones, even if among 
the conditional numbers, will be called coincident sets. This means that, for all pairs of coincident sets 16 and 
fs ), the following inequalities must be satisfied: 


fgmin< | 
min Cel. 


The number of states which are not contained in the compulsory and the conditional portions of the set will be 
called forbidden numbers, 


In using the graphical method of designing a (1,k)-pole circuit in the manner described in [1-3], the following 
operations are performed: Relay contacts are successively introduced in the circuit, the numbers are suitably con- 
verted, and the points with the coincident sets on individual circuit “sections” are combined, However, the presence 
of "direct lead - outs“ can result in alternate path circuits, which may not satisfy the assigned conditions, 


— | | 
— 

— — 
— — | 


The possibility of obtaining an alternate path can be explained by considering the network fx fa. , an) 
in a circuit of n relays Al. . An [6]. We shall assume that one of the branches in this network, S fy, at the 
point ti of section I is combined (Fig. 1) with some other network passing through the direct lead-out at the contact 
of relay Am(n=m >! 2 1) (before the contacts of relay Az have been introduced into the circuit). 


In this, besides the basic network where w = an), 
(ali. . am), and (ai. a2), also the bypass network will 
be present. 

The necessary and sufficient condition which has to be satisfied if the alternate 
. 1 path is not to form a circuit contradicting the assigned conditions is given by the fol- 

— 1 lowing inequality: 


— 
— 
7 


OY: P2P2 0 


In the particular case, this ine quality is satisfied by the te lation , 0 or 
51 ½ 0, 1. e., the presence of two unlike contacts of at least a single relay from 
Alyy,» Am in circuit or If O. the state of network vi is of no consequence for those states of the 
circuit where the network w¥,¥%,%, is closed. 


In this, if 
fx min S < | /x|, 
(5) 
network vi as well as other parallel paths from X to the input can be entirely eliminated, 


The present article offers a general number conversion method for use with the graphical method of circuit de- 
sign where bypass networks are taken into account, 


1. Operations with Sets 


For the transformations necessary in graphical circuit design, we shall introduce the following operations with 
sets that are reduced to a single base: +(addition), - (multiplication), O(combination), —(subtraction), and (inver- 
sion). If we denote by 1 the presence of a given number in the compulsory portion, by ¥, its presence in the condi- 
tional portion, and by 0 its presence in the forbidden portion, the first four of these operations are defined in Table 1, 


TABLE 1 
x v x+y xy xOy | xy" x v x+y xy xQy | x—y* 
0 0 0 0 0 6 1/2 1 4 | 1/2 4 | 1/2 
0 4/2 | 1/2 0 0 — 1 0 1 6 0 { 
0 1 1 0 0 — 1 1/2 1 1/2 1 1 
1/2 0 1/2 0 0 1/2 1 1 1 1 1 1/2 
1/2 1/2 1/2 4/2 {4/2 | 4/2 jj 


Operation x-y can be performed only if the condition |x|=|y| is satisfied. 


The inversion operation is defined by Table 2. 


Besides these operations, which are similar to three-value logic operations, we 
shall also introduce the following set transformation operations, a portion of which is 
connected with changes in the base, 


Operation P; of set division with respect to the contacts of relay Aj. By pi G). ve 
shall denote the operation of selecting from set I the numbers of those states for which 
relay Aj operates, and by BI F). we shall denote a similar operation for the case where 
the relay is inoperative, 


Thus, for instance, PsL0, 2,5, (3,6)}3=10,2/3)} 5, Py 10, 2, 5(3,6)}3=15, (%, 


TABLE 2 


x 0 | 1/2 1 


p 
2 


— 
..... 
— | 


Hence, it follows that = P,{0, 1,2,3,4,5, 6, = 12,3, 6, Thy. 
It can readily be seen that 


P. =, =0. (6) 
Operation Ht of determining the “derivatives” of numbers in passing from base Aj, . , A, to base An, 
where n>i2=1. In this operation, the number x is transformed into 2"~* numbers by adding to this number all the pos- 


sible weight combinations qj+3, di . an- For instance,for n=5, we obtain 
Hi (5, (0)) = (5, 13, 21, 29, (0, 8, 16, 24)),. 
It can be shown that 
HR (f) = Hig 
Operation Cf of reducing“ the base from Ap,. . , A1 to Aj,. . , Af. For this, the following relationship holds: 


In the case of switching over from base i to base i-1, the number of states for which relay Al operates is te- 
duced by the relay weight qj, while the remaining numbers do not change: 


for i, 


c for 1 01 (9) 


For instance, for i=3 (qj =4) 
Cz (0, 2, 4, 3, 6}, = (0, 2, 0 (3, 2)) = (0, 2, (3)}o. 
On the basis of the above definitions, it can be shown that 
Hci (10) 


3. Transformation of Sets in Circuit Design 


In the process of circuit design, each point a at the I section will be assigned three sets which will character- 
ize the requirements imposed on the network from this point to the input as well as the possibility of obtaining alter- 
nate paths, 

The “basic set“ f is a set of compulsory and conditional numbers, reduced to the base I. which indicates for 
which states the given point will be or could be connected to the input, 


The “bypass set“ & is a set of compulsory and conditional numbers that are reduced to the base n of the en- 
tire circuit, which indicate for which states the given point will be or could be connected to the input along the al- 
ternate path, 


The "forbidden set" Za is the set of the numbers reduced to the base n of the entire circuit of those states for 
which the given point must not be connected to the input. 


The following relation holds in this case: 


Cr Za]. (11) 


It follows from (11) that, for instance, the conditional portion of the basic set does not have to be assigned for 
a given forbidden set. 


For the output point x, the basic set fy corresponds to the assigned structural admittance between the given out- 
put and the common input, and the forbidden set includes all numbers that do not belong to the basic set, l. e., 
2x = while the bypass set is absent: G = 0, 


h 
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In order that the two points a and b at the section I may be combined, it is necessary and sufficient that: 
a) The sets be coincident, fa fp, and b) that the bypass set of a single point do not contain numbers belonging to 
the forbidden set of the second point, i.e., 

2.0 = O, 25 =0. 


In combining the two points a and b with the common point g (Fig. 2), the follow - 


%.%.7a 0 ing sets are assigned: the basic set 
* 
N. 26 fe = O Fo: (12) 
49% 4 the bypass set ., which is equal to the sum of the bypass sets of the points to be com- 
bined 
Fig. 2. a 


O. = + (13) 


7 the forbidden set Z,, which is equal to the sum of the forbidden sets of the points to be 
44,2, 1 combined, 


Ze = Ze — Zp. (14) 


In passing from the I th section to the (I- I)th section, the contact between relay 
A with the weight qi and the direct lead-out is added to the initial point O with the 
sets fy, Zo, and % (Fig. 3), and the sets are transformed according to the following rules, 


A. The basic sets fi and I of the opening and closing contact lead-outs are deter- 
mined and reduced to the new base 1-1: 


H= CH PI Cis Pi (fo). (15) 


B. The direct lead-out set f, is found: 


ts =f, Of, (16) 


C. In the basic sets fi and f,, the compulsory numbers belonging to the direct lead-out set are transformed into 
conditional numbers, as a result of which we obtain the following sets: 


Je, Ji = (17) 
The lead-out h, whose set will be empty, is erased in correspondence with the graphical method. 


D. The bypass sets are composed of the corresponding initial-point bypass numbers and of the new numbers 
which appear on account of networks passing through the given contact direct lead- out: 


O. = PI (®) (/s), 
= Pi (Mp) + (fs), (18) 
Dy = Dy + Ji) + (fs), 


where 
&= 1 if the lead-out with the number h remains in the circuit, 
0 if the lead-out with the number h is erased (fn is the empty set). 


E. The following forbidden sets are composed: 
Z, = P, (Ze) + P. Zo | Do), 
Z. PI (Zo) + (Zo))|Dol, 
Z. = Ly + (Zo) + 

Pi 


(19) 


— 
— — — 
— 


(18) 


(19) 


Here, the first term corresponds to the requirement that the appearance of new networks directed toward the 


circuit outputs is inadmissible, and the following terms correspond to the appearance of new bypass numbers in other 
jead- outs from the contact in question, 


Figure 4a (the number sets pertaining to each of the points xi, ti. and ui are given in Table 3) shows an exam- 


ple of synthesizing a(1,3)-pole circuit whose networks are defined by the following sets (this example was taken from 
a paper by F. Svoboda [7)): 


Ix, = (2, 3, 6}, 
IX. {0, 1. 3, 6 6, (2, 5). 
Ix, = {2, 3, 4, (ih. 


When circuits are designed by using the above method, the possibility of creating spurious bypass networks is 


eliminated, while, at each instant of time, the operations are performed only with sets of points in a single contact 
group or a single section, 


It should also be noted that, in the obtained circuits, the possibility of 
network breakdown during the transient periods of the actuation or the drop- 


A; t, ty eps A, out of individual relays is eliminated, At the same time, the method de- 
2 scribed above does not provide the possibility of using alternate paths. Neither 
4 t does this method take into account the fact that, in forming the direct lead- 

LZ u out, it may happen that not all of the possible numbers which are contained 
t Wee in the set obtained according to Eq. (16) are assigned to this lead-out; a por- 
E — tlon of these numbers can be left as compulsory numbers in sets f, and 1. 
3 5 Thus, the process of assigning numbers to the direct lead - out is not well de- 
ts fined Iii 
45 8 It should be noted that sets @ and Z cannot be composed always, but 
ty * oo only in the case where it is necessary to determine the possibility of effecting 
a combination in the presence of direct lead - outs. 
4 ty Let us now consider some additional transformations which could lead 


— to a simplified circuit, 


4. Selection of Direct Lead-Out Variants 


For selecting the set variant of the direct lead-out, we shall, in form- 
ing the latter, assign the Boolean coefficients a. .. . . to all the compul- 
sory numbers contained in its set, which is obtained by using Eq, (16). Instead 
of this, the inversions of the above coefficients (d, d:. . ) will be assigned 
to the numbers in the compulsory portions of sets 11 and . which, in corre- 
spondence with (17), must be transformed into conditional numbers, The 
Fig. 4 same coefficients will also be assigned to the corresponding numbers appear- 
7 ing in the bypass sets & and 4, and their inversions will be assigned to the 
numbers in the set 4, which is obtained according to Eqs, (18) (if a certain 
given number belongs to set db, no coefficient should be assigned to it), Similarly, coefficients 6,, .. . . are as- 
signed to the numbers for the direct lead - out at the second contact, etc. 


We shall further assign the values 0 or 1 to these coefficients, The value 1 signifies that the number to which 
a certain given coefficient is assigned remains a compulsory number, and the value 0 means that this number must be 
transferred to the conditional portion of the set (enclosed in parentheses), If the value 1 is assigned to a certain 
coefficient &; (&; indicates that the qi coefficient without inversion as well as the & coefficient with inversion can 
be taken), then the 1 value must also be assigned to all other coefficients in the section under consideration, white 
the 0 value must be assigned to the inverse coefficients, and conversely, 


The coefficient values are decided upon while comparing the sets and combinations, At the present time, there 


ate still no clearly defined rules indicating what values are to be assigned to the coefficients, and only the following 
Considerations can be given: 
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a) If one of the sets that are being compared contains such a number with a coefficient that the transfer of this 
number into the conditional portion makes the sets coincident, the 0 value must be assigned to this coefficient in 
order that the unification can be performed. 


b) If, in the case of two coincident sets, the number x in one of them is a compulsory number and in the other 
set this number has a coefficient, then the 1 value is assigned to this coefficient, which also makes unification possi- 
ble. 


o) If all the compulsory numbers of a certain lead-out have coefficients, the O value can be assigned to all 
these coefficients, which makes it possible to erase the given lead-out, 


d) If a number with a coefficient belongs to the complete set, it is convenient to assign the 1 value to this co- 
efficient. 


TABLE 3 


FES 


E 


1X1 2,3, 6 0,1, 4,5,7 t, 3, (2) 2 0,1 
X. |0,1,3,6,(2,5) 4.7 11 

X. 2, 3,4, (1) 0,5, 6.7 ue 1, (0) 2 

t 22,3 a2 0,4 us 

ts |a2 a6 {4, 5,7 — 3, (2) 7 
ty 42 42, 3, a6 0, 1,4,5, 7 8 3, (2) 7 
14 0, 511,3, 2) Ba i, Bs 2 47 0 5 
ts 2, (4) 81 5, 826 4,7 — 3, @) $7 


te 
ty 
10 


In set transformations, the assigned coefficients are also ascribed to all numbers derived from the given number, 
If, in combining, a certain number has different coefficients in the two sets, the sum of these coefficients is used 4s 
the coefficient of this number in the combined set. 


5. Use of Alternate Paths 


re This is possible if the following two conditions 
are observed in combining: 


x 2,3,8 0,4,4:5,7 2.3 2 | 0, 
X, 0, 1.3,6, (2,5) 7 * 
23,4, (i) 0,5,6,7 0,4 0,4 2 
11 | 2,3 3,(2) | 2 10,4, 
| 2 2) 0 0 3. [7 
ts 3, (2)10,1,4,5,7 1 1 6 4 
. 0, 1,3, 000.8, (1,2) 1.2 0 0 7 
ts | 2, (4) | (4,2) 0 0 5 
16 1.2 0, 3, 
t 2,3,(4)| 2,3, (1) 0 2 7 
te 0¹ 5,6,7 |_| ) 5, 
TABLE 4 
Vor 719 7 
6, (2,5) |* 
2,3, (1) J | 
0 65.8.7 
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1) Among the compulsory bypass numbers of one of the points to be combined, there is a number x (or a group 
of numbers) which is equal to the compulsory number n of output X that is connected to the second point to be com- 
bined. 


2) The given bypass number x has not been obtained as a result of the direct lead-out that is connected to the 
second of the points to be combined, 


If these conditions are observed, the number n of the second point set, obtained from the number n of the 
output X in designing the circuit, can be transformed into a conditional number if it has not been obtained by a com- 
bination with any other compulsory number, 


Thus, for establishing the possibility of using bypass networks, it is also necessary to know from which lead-out 
numbers the basic set numbers have been obtained and due to which direct lead-outs the bypass numbers have been 
obtained. 


Under the above conditions, the coefficients a, 6. . carry such information, or any indices indicating the 
appurtenance to a certain given direct lead-out can be assigned to the numbers, 


Table 4 shows sets pertinent to the circuit given in Fig, 4a, 


Since checking and comparison of sets are performed each time only for one circuit section, this method does 
not make it possible to take into account alternate paths which pass through several direct lead-outs, Thus, in the 
circuit of Fig. 4b, a network for the output Xz is formed through two direct lead-outs in the states with the numbers 2 
and 3, as is shown by arrows, The presence of such a path makes it possible to eliminate the network with the con- 


tact of relay Az. which is marked by crosses, and to leave only one opening contact of relay As between the outputs 
XI and Xz. 


6. Design of (p,k)-Pole Contact Circuits 


The graphical method can also be used for designing contact (p,k)-poles which are separable on the output side 
[1-3], whose networks are given by the (p,k)-sets fil. which correspond to networks from each output with the number 
jG=0,1,..., k 1) to each input with the number i (i=0,1,..., p=). For design purposes, the (p,k)-pole is trans- 
formed into a (1,k)-pole consisting of n+m relay contacts, where m is the smallest integer satisfying the inequality 
am > p. while each number x of the set fil is transformed into the number 


x = 42" +i (20) 


of the new set fj. The design further progresses along the usual lines until the contacts of relay Am- 1 (with the weight 
Am) are introduced in the circuit, The numbers at the obtained ends will indicate the input numbers i. 


This approach can yield positive results also in designing (1,k)-poles 
3 ＋ 7 E. under the condition that no equal compulsory numbers be contained in any 
7 a f of the output sets (the separability condition). In certain cases, a similar 

approach can be used for designing (1,2)-poles where the separability con- 
dition is absent, as is recommended in [4], For this, sets f, and fi are trans- 
formed in such a manner that each number that is present in both sets re- 
mains only in one of them. For the obtained sets fy and fi. a (1,2)-pole is 
constructed beginning from the input, for which the numbers of these sets 
are transformed into a single set according to (20), After this, the fu net- 
work is constructed, whose set comprises the above-indicated numbers as the 
compulsory numbers which are common for both sets fy and fi, while the 
numbers of the unused states and the forbidden numbers which are cot.mon 
for the f, and f, sets enter the network f set as conditional numbers, This 
network is connected between the outputs of the constructed (1,2)-pole 
(Fig. 5a). 


For instance, if the following sets are given [A]: 
ho * 0. 2, 3, 4)s, h= 1. 2. 3, 505, 
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then, by excluding the numbers 2 and 3 (this exclusion can be secured in different ways; at the present time, noclear 
indications are available about the procedure to be followed in order to obtain optimum results), we obtain 


fo (0, 2, 4)}s, 11 1. 3, 


Hence, 


for 2. 3, (6, 7005. 


Since, for the given case, m=1, we obtain the set 
7 = {0, 3, 4, 7, 8, 11},. 


The construction of this circuit is shown in Fig, 5b, the construction of network fh is shown in Fig. 5c, and the 
resulting circuit is given in Fig. 5d. 


CONCLUSION 


The graphical method makes it possible to design H-class circuits while automatically taking into account the 
unused states and bypass networks, which cannot be achieved by using the regular methods, Expetience shows that, in 
most cases, the designed circuits are simpler than the same circuits designed by other methods, while the amount of 
work involved is considerably reduced, The simplicity of operations and their identity made it possible to use this 
method for designing an original machine for synthesizing relay circuits [1,8, and 9]. 
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The present article is concerned with the algorithm and the block diagram of a specialized digital 
computer for machine tools for cutting profiles consisting of sections of straight lines and second- 
power curves, The device errors are investigated and the results obtained in accuracy investigations 
on a universal computer are given. 


The programs for metal-cutting tools with digital programming control are composed by means of universal 
computers, specialized computers, or manual calculators, 


The program obtained in these machines is transformed into unitary code, 


The extremely large volume of intermediate memory constitutes a disadvantage of programming by means of 
universal computers. Moreover, some enterprises which use machine tools with digital programming control are not 


in a position to acquire a universal computer. Therefore, a number of specialized computers for these purposes were 
developed in recent years, 


Specialized programming devices can be designed on the basis of a digital differential analyzer where binary 
multipliers are used as digital integrators, In this case, the character of the interpolation curve will be determined 
by the number of integrating units in the digital differential analyzer and by their intercoupling circuit, Examples 
of such devices are given in [1 and 3), Their disadvantage is the relatively nonuniform distribution of output pulses, 
which impairs the operation of digital follow-up systems and causes relatively large calculation errors, Thus, in re- 


producing a linear function, the error of an interpolator with a binary multiplier attains the following value under 
certain conditions: 


1 1 1 
400K 
where n 1 is the number of places in the binary multiplier [5 and 6]. 


The error of programming devices consisting of several intercoupled binary multipliers can be even larger. 


The error of programming devices where increment summators are used as integrators is much smaller. Wecan 
cite examples of programming devices with parallel increment summators, which have been developed at the Electro- 


mechanical Institute, Academy of Sciences, USSR, and sequential summators, which have been developed at the In- 
stitute of Computers, Czechoslovakia. 


Calculations of Circles, Ellipses, and Hyperbolas 
As is known, the differential equation 


z+ = 0 


can be solved by means of a digital differential analyzer where the integrators are coupled in the manner shown in 
Fig, la, 


If the pulses produced by the integrators are fed through scaling circuits to digital follow-up systems which 
move the machine tool working element along the y and x axes, this element will move along an ellipse with the 
angular velocity % = ua and the major to minor axis ratio k;/kg. Here, k1 and kg are the scaling factors of the 
scaling circuits, If these factors are so chosen that ki = kz, the working element will move along a circle. 


Similarly, the hyperbola equation %-wix=0 can be solved by means of a digital differential analyzer 
with the same circuit, but with positive feedback (Fig. ib). The working element motion trajectory is in this case a 
hyperbola whose asymptote slope is equal to ky/kg. 
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Let us find the algorithms for calculating taese curves, Consider a number of discrete points on an ellipse 


which are arranged in such a manner that the parameters of the adjacent points differ by ö. 


a b 


Fig. 1. 


Let the i point coordinates, calculated from the ellipse parametric equations, be 


2% = psin (8 4. 4), yi = © cos (0 8) 


= psin@cosd + pcosésin 4, 
71 cos 0cos 8 — £-sin@ sin ö. 
The coordinates of the adjacent (i-1)th point will be 
= psin6, 


By substituting (3) and (4) in Eqs, (1) and (2), we obtain 


11 = cos 6 -+ ayi_, sind, 


Yi = yi, cos 6 — “1 sin ö. 


141 = % cos 4+ sin ò, 
Lite = cos + sin ò, 
= cosd — 44 sin ö. 


By taking into account relation (5), we have from Eqs, (6) and (8) 


= yrcosd — ein = — sind ) cos — sind = 


= yiscos*d + ein' ö —2— sind = 2. ein d. 


Similarly, by taking into account relation (9), we obtain from Eqs, (7) and (8) 
41 = 271 sin 6. 


(1) 


(2) 


(3) 


(4) 


(5) 


(11) 


(12) 


— 
| | | 
Similarly 
(7) 
(8) 
(9) 
(10) 
— — = 
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In the same fashion, we can find the algorithm for calculating the hyperbola, 


Let the hyperbola equation be given by * A = p. We shall express the coordinates in the following 


manner: 
1 = p ch (6 + d) ch ch sh = y_,chd + ay_, shd, 
Similarly 


141 chd + ayshd, 
Lite & + sh, 
Much 0-4 sh, 
Considering that ch*§ - = 1, we obtain from Eds. (13), (14), and (16) 
= Yiu + 22 sh 6. 
From Eqs. (15), (16), and (17), we obtain 


= + Zayi4, sh &. 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


It is obvious from Eqs, (11), (12), (19), and (20) that the values of sind and ché are independent of the values 
of coordinates x and y and that they are parameters which are determined by the argument increment for each itera- 
tion cycle. By substituting w for 2siné and 2 ch &, we obtain the final definition of the algorithms for calculating the 


ellipse and hyperbola coordinates in the form of difference equations: 


Yin = — a Ti, 


Vita = + 
for the ellipse, and 


V = Yi-1 + i, 


Tits = % + 
for the hyperbola. 


(21) 


(22) 


(23) 


(24) 


For the processing of the obtained algorithms, we can suggest computer circuits with parallel and sequential in- 
crement summators. Consider a parallel circuit consisting of two registers, Rx and Ry, and two summators, Cx and Cy. 


as is shown in Fig. 2. The Rx and Ry registers consist of reversive counters, 


The y;-, value, obtained from the preceding iteration, is preserved in Ry and Cy, while the x; value, obtained 


in the preceding iteration, is preserved in R. and Ci. 


Each iteration consists of two cycles, In the first cycle, the number contained in R, is subtracted from the num- 
ber contained by Cy; during the second cycle, the number in Cx is added to the number of Ry. The ki and k units 
consist of multipliers, The multiplication factors are chosen in such a manner that kg=w/a and kia (ki and ke are 
always smaller than unity). After passing through the ki unit, the summator Cy overfilling pulses are fed to counter 
Ry and are subtracted from the number contained in the counter. After passing through the k, unit, the summator Cy 


overfilling pulses are fed to counter R, and are added to its contents, 
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Thus, iterations = Yin — — and X49 = % + are performed, The pulses leaving units ki 


and kg are the output pulses Ay and Ax, respectively, The mobile element, which is driven by the two digital follow- 
up systems that are controlled by these pulses, moves in chessboard order, l. e., it first attains the position yt; then 
Xj4g, then again yj+s, etc. In the circuit [4] constructed at the Institute of Electromechanics, Academy of Sciences, 
USSR, this movement is secured by introducing a lag line into one of the channels, 


22 

pulse 
Opening of 
fy valve A 


Shutoff of 
valve A 


Moreover, a computer circuit with sequential incre- 
ment summators can be used (Fig. 3). The circuit consists 
of two shifting registers Rx, Y and x. Ry and of other single-place summators S, which can perform the addition and 
subtraction operations, The time diagram explaining the circuit operation is shown in Fig. 4. Such a circuit was used 
in the digital computer which was developed in Czechoslovakia [2]. 


Error of the Device 


In order to determine the error magnitude, we shall first find some intermediate expressions, From (21) and 
(22), the present coordinates can be expressed in the following manner: 


5 2, 

11 20 (1 — w*) + awy4, 

Ys (1 — w*)— = (20 — ), 

4. = 0 (1 — 30 + @) + awy_, (20 - ), 


Valve A 
Shifting 
a Pulses con- 
Fig. 4. 
—— 
bz 
Fig. 2. 3 
= 2 
A 
Treaster B 
Fig. 5. 
— 
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1 
j=0 j=l 


— ꝗ ꝗ ott | = 


i 
2 ((— 1 CH) + ay. 2 1) 


This can be written as 


i 
= — 4 25 
1 = + an 


i-1 
= 
j=0 
i 
kan = 2 ((— 1)’ CH 


j=l 
From the definition, it follows that 


Zo = psing, 
cos (8 — d), 


= p sin (8 + 2i8), 
cos (0 + (2i — 1) öl. 


From Eqs, (30) and (31), we obtain 


p sin d = x,; 
cos = y_, —sind p cos = Ain ö. 
By substituting these expressions in Eqs, (32) and (33), we obtain after simplification 


Yoi_, = Y_, 008 (2i — 1) 8 — A [sind cos (2i — 1) & + sin (21 — 1) ö, 
Tai = Z, (cos 2i6 — sin & sin 210 + ay_, sin 20. 


(27) 


(28) 


(29) 


— — 
* 
(26) 
where 
: 
| | 
(31) 
and 
(33) 
——„— 
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For very small 6 values, which are actually used in practice, the latter equations can be replaced by the ap- 
proximate equations 


Yous = Y_, cos (21 — 1)8 — sin 210 (34) 
and 
Ly, = T cos (2i + 1) ö + ay_, sin 2i6. (35) 
By comparing (34) and (25) with (35) and (26), we find: 


51-1 = cos (2i — 1)ò, (36) 
= cos(2i + 1)8 
(37) 
) = sin (38) 


The Device Errors in Calculating Ellipses and Circles 


The following operation must be performed for each iteration in calculating ellipses and circles: 


@ 
= Yas 2 „ = 


In practice, however, the numbers x and y are divided into two portions, The integer portions of x and y are 
contained in N. and Ry, while the fractional portions are contained in C,—unit kz and Cy unit ky, respectively, 


During the iteration process, only the integer portions of the numbers are added, The fractions are rejected, and an 
error arises, We shall call this the rounding-off error. 


We shall denote the exact coordinate values by the capital letters X and Y, while their actual values will be 
denoted by the small letters x and V. The numbers contained in NM and Ry will be denoted by [x}jn and [y]in, while 
the numbers in Cx and Cy will be denoted by xg and Yg» respectively. Thus, we can write 


Ly + Bo» 
Yor—1 = + 


Zoi Mn + pai, 


where ¢-, and pp are errors caused by the inaccurate assignment of initial conditions, and x» and y-, are the initial 
values of the x and y variables, 


= Yas — (40) 


Xu = + df gi, 


= By (1 — + awe_, + @*2,, 4% 


Nei 


i i-1 
Mas = t an- + © 4% 2 · (43) 
The expressions for the errors (42) and (43) consist of two portions, The first two terms of these expressions rep 
resent the error arising due to the inaccuracy in assigning the initial conditions, and the last two terms are the round- 
ing-off errors. 
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In order to find the exact error value, it is necessary to know all values xof of the fractional portions for each 
iteration: 


Tog, Tzig, ig, Yous Y(2i—1) g, etc. 


This is impossible in practice, since the rejected fractional portions xi and y% will have different values for 
different initial conditions in different w; generally speaking, they can be considered as random functions, However, 
for a very large number of iterations with a sufficiently high degree of accuracy, we can use the average value of the 
rejected fractional portions instead of their instantaneous random values in error calculation, 


By denoting the average values of the fractional portions xg and yg by Xg and 78 and by substituting them in 
Eqs. (42) and (43), we obtain 


- 


Ho + ane (45) 


On the basis of (27)-(29), it can be shown that 


2 1—1 J 


By taking into account (36)-(38) and by substituting (46) and (45), we obtain 


eu = &_, cos (2i— 1) 8 — E sin 210 + ＋ sin 218 + y, [1 — cos (2i — 1) ö, 


= cos (21 + 1) ＋ ae_, sin 28 + 2, [1 — cos (2i + 1) 6) — ay, sin (48) 


If the initial values x9 and y- are given accurately, €-, and fig are equal to zero, Then, expressions (47) and 
(48) assume the following form: 


Lein + — cos (21 — 1) ö], (49) 


u, = Z, [1 — cos (21 + 1) 8] — ay, sin 218. (50) 
If we consider that ig and 7g can assume any of the possible values between zero and unity with equal proba- 


bility, 0.5 should be taken as their average value. The validity of this assumption has been confirmed by many cal- 
culations, 


In order to obtain the maximum rounding-off error, Xg = 78 = 1 (instead of 1) must be substituted in Eqs. (49) 


and (50); the maximum error value will then be twice as large as 32 calculated with respect to the arithmetical 
mean value. 


It should be noted that the maximum error actually can ne ver arise, since it is impossible to specify such con- 
ditions for which all the rejected fractional portions xg and yg are equal to unity. Thus, this estimate will definitely 


give a higher value. From an analysis of Eqs, (49) and (50) and also from the data obtained in testing a parallel cir- 
cuit, the following conclusions can be drawn: 


1, The rounding-off error is independent of the ellipse axis length or of the circle radius; it depends only on 
the value of a, l. e., the ratio of the major axis to the minor axis. 


2. For a >1, the error e is always less than 1, while the error increases with an increase in a and has a nega- 
tive value, 


) 
| 
1 
3) 
| 
nd- 
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3. For a< 1, the error j is always less than 1, while the error & increases with a decrease in a and has a posi- 
tive value, 


4. In circle calculations, the error h is always less than unity, and € attains unity for 6 “=. If, in processing 
the curve, the residues in the registers are rejected and the initial conditions are again established each time X passes 
through zero, the error will not exceed unity in processing the entire ellipse. 


5. In order to reduce the rounding-off error, it is sufficient to reduce the ig and 5g values. 


0 2n 


inte- | frac- 
R lt onal i f i f i f i f 


0 — 146 — 145 — 145 — 145 — 
00 00 74 01 00 74 74 01 00 10 
74 00 00 10 74 00 00 10 74 00 


1444 — — 
* 000 764 000 764 764 006 000 764 
764 607 000 


i 0 — | 14443 — 14442 — 
R = 5670 | z 0000 | 0000 | 5670 0013) 0000 5670 
(3000) y 5670 | 0000 | 0000 5027 5670 0000 


Remarks, 1, All numbers are given in the octuple system (decimal system numbers are in 
parentheses). 2. Here, the i value denotes the number of iterations for every 1/2. For 
instance, i= 145 for O- 36/2 to 27 for R= 60. 


For the sequential circuit, the calculations for circle computation were performed by means of the "Ural 1" 
universal computer by solving the difference equations (21) and (22) directly, The results indicate (see table) that 
the following values are obtained for 6 2 / 2 and 6 2: y=0, x=60, 100, 500, 3000, and x= O, y=60, 100, 500, 
3000, respectively (for the initial conditions x= O, y=60, 100, 500, 3000), Thus, it is not necessary to discard the 
residues in the register and to reestablish the initial conditions when y passes through zero. 


A Method for Reducing the Rounding-Off Error 


For reducing the rounding-off error, we propose an arrangement (Fig. 5) where the numbers to be carried over 
from the first places of the registers Ry and Rx are returned to the (n-) - places of registers Y and x, respectively. In 
such circuits, the rounding-off error is the smaller, the larger the scaling factor 2~!, It is obvious from Eqs. (49) and 
(50) that the rounding-off error is proportional to the value of a and to the average values Ig and Vg of the fractional 
portions, If a scaling factor larger than or equal to a is used in the proposed circuit, the average values of the frac- 
tional portions 7g and 78 will be equal to or less than 1/2a, It is obvious from Eds. (49) and (50) that p and & will 
always be less than unity, A circuit with such a correction will be called a circuit with a small rounding-off error in 
comparison with the circuit where the transfers are effected from the ith place of the Ry and Rx registers, respectively. 


Experimental Check 


The operation of the described parallel-type device was simulated in a “Ural 17 universal digital computer. 
Over ten circles with different radii and ellipses with different major and minor axes were computed, Programs were 
composed for determining the maximum positive and negative errors in calculating circles. 


The calculation data indicate that the rounding-off error for all curves can become larger than unity, This is 
explained by the fact that the author took into account the intermediate points for each iteration, Actually, the 
curves do not pass through these points due to the fact that the time interval between two successive operations in each 
iteration is extremely small, while the servosystem angular momentum is relatively large. 
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THE USE OF SEMICONDUCTOR DIODES IN MAKING 
OF SYMMETRICAL NONLINEAR RESISTANCES 


Tai Tzü-hsin 


(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 3, 
pp. 369-375, March, 1961 

Original article submitted May 27, 1960 


An analysis of the factors is made which influence the stability of symmetrical nonlinear resistances 
composed of semiconductor diodes connected in a bridge circuit, Examples are given for determina- 
tion of the sensing threshold of a magnetic amplifier using such resistances for selective rectification 
of even harmonics of output voltages, It is shown that silicon diodes are suitable for construction of 
symmetrical nonlinear resistors to be used in magnetic amplifiers of high sensitivity. 


In magnetic amplifiers in which four harmonics are rectified by means of a symmetrical nonlinear resistance 
[1,2], the latter can be made up either of carborundum resistors, or of two diodes connected in series or in parallel 
with opposite polarities, or of a number of diodes connected in a bridge scheme, The carborundum resistors are 
stable, but are subject to high losses [3], When two diodes of opposite polarities are used, it is very difficult to con- 
struct a nonlinear resistance having a good symmetry of the voltampere characteristic around a selected point, be- 
cause the individual diode specimens differ a great deal in their characteristics, The symmetrical nonlinear resist- 
ance produced by four electronic valves connected to a common cutoff voltage (Fig. 1) [4] also has its shortcomings; 
these are: first, the necessity of employing for the cutoff a battery or other source of constant de voltage; second, the 
great losses introduced by the high resistance of this voltage source, especially if it is of low power; and third, the 
eventual presence of a so-called “zero current,” that is,a current in the output when no voltage is applied to the a, 
b terminals of such a symmetrical nonlinear resistance. 


Of great interest is a bridge scheme (Fig. 2) which has in its diagonal branch a 
“stabilitron” voltage regulator consisting of a silicon breakdown or zener diode of one 
of the D808-D813 types. Potentiometer rheostat or slide wire Rp serves to adjust the 
b bridge to symmetrical output. Such a bridge has low impedance and small physical di- 
mensions, In obtaining of symmetrical nonlinear resistances of high stability, the deter- 
mining factors are the magnitude and stability of the forward resistances and of the re · 
Fig. 1. verse currents (or reverse resistances) of the rectifying diodes in the bridge arms. 


“It has been proven that at higher voltages—say, above 5-6 v— an avalanche breakdown rather than the Zener phe- 

nomenon is responsible for the operation; therefore, there is at present among semiconductor engineers a tendency to 
use the term “breakdown diode,” rather than “zener diode,” as the generic term. In this translation, as well as in the 
original Russian text, the name “zener diode” is not used [Translator’s note]. 
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In Tables 1 and 2 are presented characteristics of germanium diodes of type DGTs-27 and of silicon diodes of 
type D202, all obtained experimentally, The diodes were chosen to have peak inverse voltages ujny of 16 v. which 
is much higher than the highest possible breakdown voltage in the bridge diagonal; in that case the reverse voltages 
normally applied to each of the rectifying diodes in the bridge arms are approximately equal to the rated voltage out- 
put from the type D808-813 silicon st abi litron. The data of Table 1 show that the reverse currents with a potential 
of 16v applied to DGTs-27 germanium diodes reach tens and hundreds of microamperes at elevated temperatures, 
Farther below in this paper it will be shown that in a typical magnetic amplifier for measurement purposes which is 
to have sufficient sensitivity for sensing an input of 10-“w, the highest permissible "zero" current (output current 
with zero input to the amplifier) must be less than 3.5 ya, The high reverse currents and the large scattering of the 
characteristics of DGTs-27 germanium diodes make it impossible to use them for construction of a bridge of the type 
under discussion. In addition, in the DGTs-27 diodes after application of the reverse potential, a change of the re- 
verse current with time is observed as illustrated in Fig. 3. This phenomenon is especially conspicuous at elevated 
temperatures, At 55°C also a lack of constancy, or reproducibility of the reverse currents in the very same DGTs-27 
diode in successive tests, is observed: the reverse current varies by several microamperes from test to test. 


TABLE 1 
Characteristics of DGTs-27 Germanium Diodes 
d re- 
verse current with| Forward verse Current w Forward — 
No. of |16v applied reverse sistance fo No. of |16v applied reverse g ma form ard 
diode voltage, pa current ohms] diode voltage, pa current ohms 
20°C 55°C 20°C 20°C 55°C 20°C 
1 8.9 92 92.4 6 17.8 134 92.0 
2 14.8 109 99.0 7 18.0 240 97.3 
3 20.0 256 90.5 8 19.4 126 91.3 
4 15.2 193 101.0 9 14,8 134 98.2 
5 10.0 74 94.1 10 12.0 184 93.0 


TABLE 2 
Characteristics of D202 Silicon Diodes 


verse current with 16v | Forward resistance (ohms) at forward current of 
pplied in reverse direc- 
* tion, pa 2ma | 20 ma 

2˙ [2050 20°C 90°C 
1 0.05 0.29 322 264.5 37.8 33.55 
2 0.00 0.40 318 260.0 36.8 32.45 
3 0.00 0.41 320 261.0 37.1 32.85 
4 0.05 0.70 320 270.5 38.4 34.40 
5 0.00 1.85 310 251.0 36.2 36.90 
6 0.00 0.15 318 259.5 37.1 33.25 
7 0.20 2.88 318 262.5 37.0 33.40 

8 0.10 1.90 320 37.5 


Let us consider the possibility of using silicon diodes of type D202 in this bridge. The forward resistance of 
these diodes is about three times as large as that of the DGTs-27 germanium diodes, but it varies much less from 
specimen to specimen at a given temperature. The reverse current in the D202 diodes is very small even at 90°C 
(Table 2); however, here the forward resistance of the diodes is of much more importance than the reverse current: 
it is the forward resistance which practically determines the sensing threshold (minimum of detectability) in a mag- 


netic amplifier when the bridge discussed here is used. In view of the above, ultimately silicon diodes of types D202- 
205 were chosen as rectifiers for the arms of our bridge. 


| 

12 


202- 


In Fig. 4 a static voltampere curve of the bridge is shown when the silicon diodes were used in the bridge arms, 
and a silicon “stabilitron” in the bridge diagonal, The dynamic resistance at a current of 20 ma is 
Au 11.5 — 10,6 ohm. : 
R= (20 —0.11) 10° 45.1 
The breakdown diode in the bridge diagonal has no influence on the symmetry of the bridge when all four 
diodes in the bridge arms are equal, because the bridge is being opened for passage of current at the same applied 


voltage regardless of the latter's polarity, Sources of asymmetry are inequalities of forward resistances and of re- 
verse currents of the rectifying diodes in the bridge arms, 


5 Muy 
5+ 
Fig. 4. 
u 
1 TL: 
Te 
Ha 1 
Fig. 5. 
7 4 Let us first examine how the bridge symmetry is af- 
Ur %, fected by the forward resistances of the arms, or, in other 
words, how the sensing threshold of a magnetic amplifier 
50 using such a bridge is affected by these arm resistances, 
Fig. 3 We assume that a rectangular ac voltage as shown in 


Fig. 5 is applied to the a ande terminals of the bridge cir- 
cuit of Fig. 2. The character of operations of the circuit of 


Fig. 2 is not changed thereby, 
During the positive halfperiod the total current passed through the bridge (see Fig. 2) is 
pos it Tin + (1) 


where lif and ig¢ are the forward currents through diodes 1 and 3, and dur and i - the reverse currents through diodes 
4 and 2, 


It is assumed that during the halfperiod under consideration the currents in diagonally opposite arms of the 
bridge are equal, that is, that 


ig lar (2) 
and 
— (uy + Ry) 


where u is the voltage applied to the bridge during the period of time t between 0 and T/2(0=t=T/2), uo ~ the 
breakdown blocking voltage of the diode in the bridge diagonal, Ryf and Rgf— the forward resistances of diodes 1 and 3, 
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R the sum of all active resistances in the output circuit including Rp pos. which is the resistance of the part of po- 
tentiometer Rp in series with Ref still inside the bridge (Fig. 2). 


When 1,Rq << u- up, equation (3) is simplified to the following: 


— 


ig (4) 


From (4), (2), and (1) we obtain the total current through the bridge at OS ts T/2: 


2 


In a similar manner the total current through the bridge during the negative halſperiod (IT/ 2 5 t= I) is obtained 
which is 


= ix = iu + 


Ineg = tx in +73 (6) 


where iir and iar are the reverse currents through diodes 1 and 3, and ig¢ and igs - the forward currents through diodes 
2 and 4, 


In order that the bridge be symmetrical, the following condition must be fulfilled: 


Ry ＋ RH = + Re. (7) 


When the condition of equation (7) is fulfilled, and in addition the reverse currents in the arms are equal, the 
de component for the entire cycle is equal to zero, that is, 


u—t, u — 


0 
=la A FR — (u + 71) — 


Let us now assume that the forward resistances of three diodes are equal, but the resistance of arm one has 
changed by a quantity ARif. 


When there is a change in the forward resistance, the reverse voltages applied to the diodes are changed only 
insignificantly; therefore, it is possible to disregard the changes in the reverse currents caused by the resistance 
change ARIf. Then, for the deviation of the “zero” current from actual zero, that is, ſor the dc component, we have 


u — u,, (u —u,) ANif (8) 


By means of (8) it is possible to find the highest permissible variation in the forward resistances of the bridge 
arm diodes corresponding to any given acceptable sensitivity level. 


When to an arrangement according to the circuit diagram of Fig, 6 a symmetrical sinusoidal potential is ap- 
plied which is sufficiently High for magnetic saturation of the cores during the major part of each halfperiod, then 
the output voltages have the shapes of pulses [I]. For the sake of simplicity, let us assume that these pulses have 
rectangular shapes, that is, that the potentials at the terminals of the output winding maintain constant values for the 
entire duration of the halfperiod Tp of the pulse (Fig. 7). In this case quantity Aly is determined by the following 
formula: 


(u—u,) ARs 
+ Ry + (9) 
From this it follows that, to have a desired threshold Alo of current sensitivity, the permissible variation of the 
forward resistance is 
2 T (Ry + — 
ARy =Al, 0 N * (10) 
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Let us now examine how the reverse currents in the bridge arms affect the symmetry of the bridge. We make 
the following assumptions; 1) the reverse currents of three diodes remain unchanged and are all equal to each other, 
while the reverse current of the fourth diode is changed by Alu: 2) the change of the forward resistances produced by 
the change Aigr can be disregarded, Throughout T/ 2 ts T the total Ineg current through the bridge remains un- 
changed at the value specified by equation (6), The change of the reverse current of the fourth diode which takes 
place at OS ts T/2 causes a deviation of the zero,“ that ls, the appearance of a constant de component Aly in the 
output of the bridge. Because, when ig is changed, i;f is also changed, the zero deviation is 


Al, = Ipos —Ipos = (ig +i )—(ia +if ) = Aig — Aig. (11) 


Let us find the increment Aiy¢ of the forward current of diode one which is caused by the change of the reverse 
current lar 


C 
Fig. 6. Fig. 7. 


When the reverse currents are changed, the voltage at the bridge terminals a, b remains practically unchanged, 
With constant voltage applied to the bridge terminals a- b, the following changes take place simultaneously: 
1, The reverse current through diode 4 increases from ig, to ig; the voltage applied to it decreases by Aug, 

2. The forward voltage across diode 3 increases by Aug, the forward current through it by Aig. 


3. The forward voltage across diode one is decreased by Au; and, assuming that the voltage across the diode in 
the bridge diagonal remains constant and equal to uo, the forward current through diode one is decreased by Alf 
= tif. 


4. The reverse voltage across diode 2 is increased by Aug, but the current through it remains unchanged, 
The reverse voltage across the diodes is in a section of the characteristics where the reverse current is inde- 


pendent of the applied reverse voltage. Therefore, when the reverse voltage applied to diode 2 is increased by Aug, 
the reverse current remains unchanged, 


The current through the diode in the bridge diagonal is is decreased by an amount Ai, which is equal to the 
change of the current in diode one (according to the first networks law of Kirchhoff applied to the currents meeting 
at point c, Fig. 2). For point d (Fig. 2) we have 


Ai |=|Aig |+ | |. (12a) 


When the voltage across the a- b terminals of the bridge is constant, and the breakdown voltage uo of the 
diode in the bridge diagonal also is constant, the voltage increase across diode one is equal in absolute value to the 
voltage decrease across diode 3, that is, [Au] = Au. Therefore, since diodes one and 3 have identical character- 
istics, the variations of their forward currents are also equal in absolute value, that is, 


| =| Adige |. (12b) 


From (12a) and (12b) we have 
[Aig = 2 Ai |. (12¢) 


Substituting (12c) into (11), we obtain the following expression for the zero shift caused by a change in the re- 
verse current of one of the four bridge arms diodes: 


Ai 
Al, = (13) 


Thus, we see that the change in the total bridge current is smaller than the change of the reverse current of one 
of the diodes. 


Temperature changes do not cause an asymmetry of the bridge, provided the following equation prevails: 


dr — tar )+ d — ier) = (i — ier) ＋ (ir — iir). (14) 
wherein the i's marked with apostrophes designate the currents after the temperature change. 


For greater convenience in selection of the diodes, equation (14) may be split up into the following two equa- 
tions: 


ig in = iat + ty (15) 
ig Tir = in + ip. 
(16) 


It is recommended to select the diodes to fit equation (16) for an elevated temperature, since in that case usual- 
ly equation (15) too is fulfilled as a natural consequence, 


In conclusion it can be stated that with respect to the reverse currents the diodes must be selected in such a 
way that, first, the sum of the reverse currents in one pair of diametrically opposite bridge arms is equal to the sum 
of the reverse currents in the other pair of bridge arms, and, second, that the highest values of the reverse currents oc- 
curring within the range of operational temperatures are taken into account for establishment of this equality. 


TABLE 3 


lee, ohms 2 ols ohms ohms 


400 | 200 | 200 


Let us evolve the relationships between the sensing threshold and the "zero" current when the forward resistance 
in one of the arms of a nonlinear symmetrical bridge, connected to the output of a magnetic amplifier (Fig. 6), is 
changed. 


Designating the lower threshold of sensing (minimum detectable quantity) of the magnetic amplifier, in terms 
of power, by P, and the shift in “zero” current at the amplifier output by Alo, we have [I] 


where ky is the current amplification and Re- the sum of all active resistances in the control circuit. 


From (10) and (17) we obtain an expression for the highest admissible variation of the forward resistance for a 
given sensing threshold, which is as follows: 


JJ. 
| ky | late wa ARs, ohms 
wo | | | 20 | 70 | | 0.38 | 0.35 
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From (17) and (18) we can find for any required sensing threshold the corresponding highest permissible varia- 


tion in “zero” current AI, and the corresponding highest permissible variation of the forward resistance. For example, when 
typical specification data for an instrumental magnetic amplifier with low sensing thresholds, built according to the 


diagram of Fig. 6, are substituted into formulas (17) and (18), the data of Table 3 are obtained, It was assumed that 


e- Ty= 1/20, u- 20 v. and uo 10 v. 
From the calculated values of permissible increments of the zero“ current given in Table 3 and the data pre · 
sented in Tables and 2 we come to the conclusion that using germanium diodes it is impossible to create a sym- 
13) metrical nonlinear resistance whose zero“ current would be smaller than 3.5 pa, 
- CONCLUSIONS 
It is shown that inclusion of a stabilitron voltage regulator, consisting of a breakdown silicon diode, in the 
bridge diagonal makes it possible to obtain nonlinear resistances having uneven voltampere characteristics, The cus- 
1 tomary use of germanium diodes in the bridge arms is very little suitable for creation of symmetrical nonlinear re- 
14) sistances of high stability, because of the high and variable reverse currents of these diodes, especially at elevated 
temperatures, Silicon diodes of D202-205 types have high reverse resistances even at 90°C, The forward resistance 
of these diodes is about three times higher than those of the germanium diodes, but it is much more uniform and 
F invariable at a fixed temperature, Their use makes it possible to build magnetic amplifiers which.operate on the 
principle of rectification of even harmonics, by means of symmetrical nonlinear rectifiers which possess low sensing 
(15) thresholds at temperatures varying within wide limits. 
The author expresses his thanks to M, A, Rozenblat for a number of valuable suggestions and comments, 
(16) 
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SENSITIVE MAGNETIC MODULATOR WITH TWO-PHASE INPUT 
F. I. Kerbnikov and M. A. Rozenblat 
(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 3, 
pp. 376-382, March, 1961 
tance Original article submitted June 1, 1960 
). is 
Circuit and theory are presented of a two-phase magnetic modulator which uses perpendicular to 
= each other magnetic fields and is supplied with power from a source of rectangular voltage pulses, 
It is shown that such a modulator has substantial advantages over other, more widely used types of 
modulators. 
(17) Results of an experimental investigation of this modulator also are presented, 
It is known that of all types of magnetic modulators the modulators operating on the principle of frequency 
* doubling have the lowest sensing threshold [1]. The limit of sensitivity of such modulators is set by the magnetic 
noises and is within the range of 10°” to 10°" w, 
(18) However, such a low sensing threshold can be assured only by the use of complicated schemes for removal of 


even harmonics from the supply voltage and suppression of the disturbance voltages (uneven harmonics) in the modu- 


— .. 
20) 
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lator output, Another shortcoming of the frequency doubling modulator is the necessity of connecting into the circuit 
a blocking filter or a large inductance to eliminate the shunting action of the low resistance signal source with re- 
spect to the modulator output. 


Use of modulators with perpendicular to each other fields [2] makes it possible to some extent to eliminate the 


mentioned shortcomings of modulators with frequency doubling by substantially lowering the parasitic disturbance 
voltages in the modulator output, 


There are further improvements possible in the characteristics of magnetic modulators using frequency doubling 
such as; decrease of the disturbance voltages, lowering of the sensing threshold, improvement of the curve shape of 
the output voltage, and elimination of the necessity of having a blocking filter in the control circuit loop; these im- 
provements can be achieved by the use of two similar modulators each with two perpendicular to each other fields, 
supplied by a two-phase voltage of rectangular shape, 


The present article is devoted to a presentation of the results of an investigation of the modulation scheme using 
a two-phase supply voltage, suggested by the authors, 


1. Scheme of the Modulator 


The magnetic modulator consists of two ferrite cores 1 (Fig. 1), each of which is provided with an excitation 
winding W, which serves for production of the transverse alternating magnetic field, a control winding We, and an 
output winding W. [2). 


As an electric power supply, a two-phase magnetic transistor generator-converter is used, consisting of two 
identical de to ac converters and an arrangement for synchronization of these two converters [3]. 


This scheme (Fig, 2) is distinguished by the feature that there are in the output two rectangularly shaped voltages 
U4, Us which differ from each other by 90° in phase, The generator has high efficiency, and its output voltage has only 


a very small content of even harmonics, Moreover, these harmonics can be practically reduced to zero by adjustment 
of a variable resistor Ro. 


7 


W, 


The rectangular shape of the voltage supply pulses makes it possible, as will be shown, to obtain practically 
sinusoidally shaped voltages in the modulator output, 


In the scheme of Fig. 2 the voltages of doubled frequency, which are produced by induction in the control wind- 
ings We when a signal is present, mutually compensate each other. Therefore, there is no need any more for a spe- 
cial filter or inductance in the control circuit loop for elimination of these doubie frequency voltages. 


2. Operation of the Modulator When in the Idling State 


Let us examine one of the excitation circuits of the modulator shown in Fig. 2, when to the excitation winding 
Wi an excitation potential of rectangular shape is applied, while in the control winding W, there is only a very small 
signal applied, 


During one halfperiod of the supply voltage (0<t<T/2), the process in the excitation circuit is described by 
the equation 


W,S,10* + Ri, =U, 0 


— 

[wi 7 
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where U, Is the amplitude of the supply voltage, R- the resistance of the WI winding, Sy - the cross-sectional area for 
the excitation flux, 


When the resistance of the excitation circuit loop is sufficiently small, the i component in equation (1) can 
be considered negligible for the major part of the halfperiod under examination, Neglecting Nui, integration of (1) 
gives 


=—By+ (2) 


wherein - By is the value of the magnetic induction By at the moment t=0, 


iges 
ent 
Fig. 2. 
It follows from the symmetry condition that, when a steady state of operation has been reached, at t= T / 2 there 
must be By B. Hence we have 
B. 10 10 U 
11 4 (3) 
where f is the frequency of the excitation voltage. 
Considering the material of the core to be magnetically isotropic we have 
H, HI I 
nd- where B= VBI + Bi, H = Hi + Hi. By and Hy here are respectively the induction and field strength 
* (magnetizing force) in the longitudinal direction, and By, and Hy - the induction and field intensity in the transverse 
direction, 
Under the influence of a signal ig in the control winding, a longitudinal field strength Hg= O. 4r iow, / I is pro- 
8 duced, to which an induction in the longitudinal direction corresponds which is 
all 


B 
By H. Hop, (4) 
where p is the static magnetic permeability of the core material. 
Under these conditions on the terminal of output winding of the modulator under examination the following 
emf will appear: 


= — 10 H, (5) 


˖˖* 

— 

| 


We must point out that the n used here must be determined from measurements on the so-called “average” mag- 
netization curve, which is the curve located midway between the ascending and descending branches of the hysteresis 


loop. 


If we limit our examination only to small signals, we can assume that the law according to which p varies with 
time is univalently determined exclusively by the variation with time of the transverse inductance component BI, and 
practically does not depend at all upon the magnitude of the signal. Therefore, expression (5) can be rewritten as fol- 


lows: 


e, = —W,S,H, 10% (6) 


where B=B,. 


aB dt’ 


If we neglect the potential drop across the ohmic resistance of the excitation circuit, then for By formula (2) is 


valid, and from it we obtain 


dB 10 
a (7) 
Substituting this into (6) we have. 


Let us assume that the average magnetization curve of the core has been obtained in graphic form by means of 
suitable measurements (Fig. 3a). It is easy to plot from measurements made on this curve the curves expressing the 
quantities = / H. = f G) and dy/dB =) as functions of B (Fig. 3b). 
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If the condition corresponding to equations (2) and (8) prevails, 
that is, when the ohmic resistance of the excitation loop is practical- 
ly zero, the curve of the output of the modulator ez will coincide with 
the curve dp /B =) within the limits from B= HY to B = = 
Experience has shown that it is not difficult, when using ferrite cores, 
to select such a value for By (by a suitable choice of supply voltage), 
that ) dy/dB will be close to a harmonic function, 


That means that within the indicated range of B, the relation- 
ship between B and h can practically be represented as follows: 
— Emin — "min (9) 
1 
where B is determined by formula (2), B% by formula (3), and n min 


is the value of the magnetic permeability at B= A5 Substituting 
this value of u into (8) and taking into account (3), we obtain 


= — H, (Umax — sin 20t. (10) 
Thus, for determination of the amplitude of the output voltage 


it is sufficient to know the maximum value of the permeability from 
the average magnetization curve, and the minimum value of permea- 


‘bility u min which corresponds to B= 45 Since usually u min max: 


for approximate calculations it is sufficient to know the value of max · 


Let us point out that the presence of a finite ohmic resistance 
in the excitation circuit has an effect only when B approaches the 
values of 45, that is, in the region of saturation of the core. More- 


over, the ohmic resistance plays a part favoring a sinusoidal shape for the output voltage ez, since because of it at 
B= 45% practically dB/dt * o. Therefore, even in the case that the derivative du / dBß + O at B= 4Bo, the emf ezbe ; 
comes equal to 0 because dB/dt o [see formula (6)). 
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In the next halfperiod when the supply voltage changes from +U,to Ui. we obtain for the emf in the output of 
the first modulator core an expression which is exactly the same as expression (10) obtained for the first halfperiod, 
because p (t+ T/2) = p(t). 


Taking into account the 2/2 phase shift between the supply voltages of the two core windings, we obtain the 
follow ing expression for the output voltage from the second core of the modulator: 


es (wt) = — (ot + 2 o/ Hie (uman Heir) sin 20f. (11) 


When the output windings of the two modulator cores are connected with each other in series opposition the 
overall output voltage is 


fut = = — Hy 10° —Pmin) sin (12) 
Since Uc = Nic we have for the voltage amplification the following expression: 
Your out 


de ie Re" 
where Uout is the amplitude of the output voltage. 


Substituting Ho 0.47 ic / I into (12) and (13), we obtain for the amplitude of the modulator output voltage, 
and for the voltage amplification: 


3. Hees,, umax — hmin 
Uour = 17 10 2 
and 
3. 22. SI Bmax — Mmin 
Ky = ~ Rely 10° 2 (15) 


where I is the average length of the core. 
3. Time Constant of the Modulator 


Under the influence of a signal in the core a dc component of the induction is established which is. 


Qn Qn 
1 
Boy = => vl. au. (16) 


Substituting p from (9) and carrying out the integration, we find that By parallels gH», where pig is the mag- 
netic permeability of the core for the de induction component, which is determined by the formula 


py = Pin, 


The inductance of the control winding is 


0.8 S — 
Le = — mia. (18) 
The time constant of the control circuit can be calculated by means of the fundamental formula for charge or 
discharge of an inductance: T Læ Re where Re is the overall resistance of the control circuit, 


4. Operation of the Modulator under Load 


An exact account of the effects of the load resistance leads to cumbersome calculations, However, for the type 
of modulator under discussion an approximate account of its behavior under a load R is possible by the use of the 
follow ing simplified method [4], The modulator is considered to be a generator of anemf=Upgy of doubled fre- 
quency having an internal inductive impedance Z;=2wL;. In this approximate computation we neglect the variation 
Of Uput and Li with magnitude of load current. The LI considered here is the average value during a period of the 
inductance of the output winding Wg, and therefore can be calculated in a manner similar to that of the inductance 
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0.8 Umax Emin 
ly ort 2 


The amplitude of the second harmonic in the voltage across the load resistance is 


Vout 
U = 
V + (20L,)* 


Taking into account (14), (18), and (20) we have 
1.6 Sy Umax + Mid 


1 105 


It is evident that the output power reaches its maximum when the load resistance is 


5. Results of Experimental Investigation 


An experimental investigation of a modulator with cores as described was carried out for the purpose of com- 
paring the actual characteristics with the theoretical, and for determination of the lowest sensing threshold. The re- 
sults of this investigation are presented in the accompanying table. 


— B., gauss Sn. 


6700 1100 2670 0.68 


The windings on the cores of the modulator which was used in these 
tests had the following parameters: W;= 100, WZ 1500, We =1000 turns; N- 
= 155 ohms, 


The oscillograms of the output voltage are shown in Fig. 4. The oscil- 
logram in Fig. 4a is the one obtained from a single core with a signal of 
WW 10%, and the oscillogram in Fig, 4b is the one obtained when both cores 
were used, It can be seen that the shape of the output voltage is already prac- 
tically sinusoidal when a single-phase modulator is used; with a two-phase 
modulator, the amplitude of the output voltage is twice as large, and -_ con- 
tent of higher harmonics is decreased. 


The voltage amplification Ku. computed by formula (15) for the idling 
state, was 386. The maximum value of voltage amplification Ky actually 
measured at f =500cps with a high resistance load, in idling condition, was 368. The ver? of the parasitic second 
harmonic in the modulator output was found to equal an input signal of the order of 4-107 


As is well known [2], a single-phase modulator with transverse fields makes it possible to obtain a sensing thresh 
old of the order of 10 “w. In Fig. 5 are presented the experimental curve (dashed line) and the curve computed by 
means of formula (21) (drawn-out line) for the relationship between the load resistance and the voltage across it at à 
current of 1· 10˙ amp on the control winding. The experimental curve practically coincides with the theoretical 
curve, 


Fig. 4. 
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CONCLUSIONS 


The proposed scheme of a modulator with two transverse fields using two rectangular pulses shifted in phase by 
90° for voltage supply makes it possible: a) to obtain lower thresholds of detectable power, that is,higher sensitivity, 
than by a modulator with a single core; b) considerably to simplify the circuit of the amplifier in which the magnet - 
ic modulator is used, since there is no more need for filters in the control and power supply circuits; and c) to obtain 
a practically pure sinusoidal shape of the output voltage. 
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CONTACTLESS SWITCHES WITH STEEL £310 MAGNETIC CORES 


K. G. Mityushkin 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No, 3, 
pp. 383-392, March, 1961 

Original article submitted May 12, 1960 


Magnetic contactless pulse switches (distributors) based on standard Steel £310 (KhVP) cores are de- 
scribed, The calculated characteristics and parameters of individual single- and double-throw dis- 
tributors are given. 


Pulse switches (distributors) based on square - hyste resis - loop (S Hl.) ferromagnetic components are presently used 
to an ever increasing extent in contactless automation and remote control devices [1-4], Either metallic daminated) 
cores made of special alloys (types N50P and N65P) or SHL ferrite cores are used in these components, The magneti 
characteristics of SHL materials that are most widely used are given in Table 1. 


TABLE 1 


Basic characteristics 
By, gauss Kp =B,/B, 


Ferrite 1450 0. 96 
2000 0,92 
2500 0,93 
2300 0.93 
2700 0.93 

15000 0. 98 

13500 0. 98 


Core type 


This article is concerned with the possibility of equipping contactless switches with standard laminated toroidal 
cores that ure made of Steel £310 (GOST 5021-55), which are widely used in various electronic devices as intermedi- 
ate and output low-frequency transformers. 


The production technology of these cores is much simpler than that required for SHL alloys, and their cost is 
considerably lower. 


The hysteresis loop of a Steel £310 core which was plotted for 50-cps alternating current is shown in Fig. 1. 
Although the £310 alloy does not belong to the group of SHL alloys, experiments indicate that these cores can be 
used in all known pulse distributor circuits which use SHL magnetic components, In this, the following advantages of 
distributors with laminated cores are preserved: the high specific power of the output pulse per unit volume of the 
material, the high ternperature stability, and the possibility of making cores of any size required in electrocommuni- 
cations practice (the outside diameter of anchor rings can attain a few hundreds of millimeters), Due to this, it is 
possible to obtain pulses with a considerable amount of energy at the distributor (switch) output, which often makes 
it unnecessary to use the intermediate amplifying units that connect the distributor to the load. 


Due to the nonrectangularity of the hysteresis loop of £310 cores, compensating elements, which serve forcom- 
pensating the noise arising in the working components, must be used in distributor circuits. 


n should be noted that the maximum outside diameter of SHL ferrite rings that are produced by the industry does 
not exceed 10mm. For technological reasons, this diameter cannot be greatly increased, 

**Emf noise arises due to departures from the rectangular shape of the hysteresis characteristic in shifting the operat 
ing point along saturation segments 0-0" and 1-1" (Fig. 1). 


ö²1—. . 
Core 
— designation Ho, oe 

1.1 
0.3 
0.8 
1,25 
0.15 
0.1 

0.05 


edi- 


The circuit of the compensating element CE is shown in Fig. 2. The compensating element does not have a 
control winding, and, for each operating cycle of the circuit, it produces noise pulses which are in opposition to the 
noise pulses of the working element WE. The same figure shows an oscillogram of pulses in the load circuit in the 
absence (a) and presence (b) of the compensating element. 


Common compensating elements are used for noise compensation in pulse distributor (PD) circuits, 


Figure 3 shows a double-throw PD circuit with two compensating cores, where one core is provided for each 
phase. The compensating windings are connected to the coupling circuits between the elements and to the load cir- 


cuits, Element CE, compensates the noise from all odd WE elements, while element CE, compensates the noise from 
all even WE elements, 


The single-phase PD circuit (Fig. 4) requires a single compensating element, whose 
10° B.gauss windings are connected to the charging circuits of the coupling capacitors and to the load 
circuits, 


The engineering calculation of the parameters of a single magnetic element withan 
£310 core can be performed under the assumption that the resulting hysteresis loop of the 
[ working core is strictly rectangular if compensation is provided (thin line contour in Fig. D. 


3 Such an assumption provides a very simple approach to the solution of this problem, 
16 26h, and it yields results that are satisfactory for practical purposes. 


For the circuit of Fig. 5, the following equations hold under the assumption that the 
core magnetization reversal follows the path from state 1 to state 0 under the action of a 
* square pulse of voltage U: 


dB 
Oy = + 10%, (1) 
Fig. 1. 


el 10 =i, N. (2) 


— Hel. (3) 
The initial core state is characterized by the point 1, and its final state is given by point 0 (Fig. 5). 


By simultaneously solving this system of equations, we obtain the following expression for the variation of mag- 
netic induction in time: 


B. = Hel (4 — 


10-8 (: + (4) 


where K- ale isthe load resistance, which is reduced to the working cireuit,kjg K is the forcing co- 
L 
efficient under the operating conditions, and I, = — ° 


w 
The total core magnetization reversal time (from point 1 to the point 0) is found from (1) for By = Br. 


2B, 10-* (1 + 
— 1) 


It is clear that, for a total magnetization reversal in the core, it is necessary that the operating pulse time tp 
be not less than t,. 


According to (2), the load current amplitude is equal to 


(5) 
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Fig. 2. W,,) working windings; Woo) control winding; WI and 


We) load and compensating winding, respectively; R) load re- 
sistance, 


b 
— 
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— 
— —ü— 
From pulse generator 
Fig. 3. 
| 
From pulse generator 
Fig, 4. 
Uy 
Fig. 5. 


By substituting in this expression the value of dB/dt, which is found from (4), we obtain the following expres- 
sion after some transformations: 


Hel 


(1 (7) 


The current through the ae winding during the magnetization reversal process is found from (3) by taking 
into account (7): 


= INI. +1]. 
1 . (8) 
The working circuit current consists of two components: the no-load current Il and the load current IL, which 


is reduced to the working circuit, After the magnetization reversal is completed, the working circuit current becomes 
equal to 


Uy 
4. (9) 


The diagrams of the magnetic induction and of the load and working circuit currents are given in Fig. 6. 


The satisfactory agreement between experimental and calculated data for a single magnetic element with 


compensation makes it possible to usethe well-known method for calculating the parameters of pulse distributors for 
SHL magnetic elements [3 and 4) when £310 cores are used in these circuits, 


Example 1, By using the above equations, we shall determine the parameters of a magnetic element that se- 
cures a current equal to I; = 1 amp in the load winding for the resistance Ry =5 ohms and a pulse duration of 5 msec, 


The core dimensions are 30 x 20x 6mm, the remanent induction is Br= 12,000 gauss, and He = 0.8 amp/em. 
The core cross-sectional area is $=0,30cm*, and the median line length is 2 =7,8cm. 


By substituting in Eq. (5) the expression for the difference kfg —1 found from (7), we obtain 


25, 118 


t, J] LAL 


Hence, WI = 350. 
We shall find the forced-conditions coefficient from expression (7): 


(1428) 


* = 


For a pulse source that is close to a current source, Ry RL. Therefore, the expression in brackets can be as- 
sumed to be equal to unity. Then, kfo * 56. 


The operating current amplitude is found from expression (8), By taking into account that Rj, & Ry, we have 


Hd 
ly = 3.5 amp. 


In calculating ly, the number of turns in the operating winding was assumed to be 100, 


As was shown by suitable experiments, the pulse in the load winding of an element with an £310 core of this 
size is equal to 0.95 amp for wy =350 and wy =100 and a duration of 5,1 msec, 
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Example 2. By using the method described in [3], we shall determine the parameters of a double-throw pulse 
distributor that is loaded directly with an output RES-6 electromagnetic relay, The schematic diagram is shown in 


Fig. 7. 
The output relay data are the following: The actuation current is 70 ma, the winding resistance is R= 135 ohms, 
and the actuation time for the voltage U; = 18 v is 7.8msec. The actuation factor of assurance for Uy, =18v is k,. = 
= 2,05, 
Consequently, the output pulse from the distributor ele- 
ment must have a power of the order of 2.5w and a durationof 
3 approximate ly 8. O msec, 


For determining the core dimensions, we shall find the 
product wis from expression (10): 


t 
wus = = 600. 


4 — te — We shall choose two 30x 20 X 6mm cores with an over- 


all cross-sectional area of S=0.6cm*, Then, wi = 1000 turns, 


i) T-- + The power Py speat on the core magnetization reversal 
hi during its actuation is calculated by using the equation 


Iwm 


„ 


(11) 


x” 


Fig. 6. U) Operating pulse source voltage; Bt) core 
induction; i,) load circuit current; i) working cir- where V is the magnetic core volume (cm*), and A, is the 


cuit current, hysteresis loop area (gauss =P). The values of V and Apare 


equal to 
V= Slay = 4.7 cm’, 


Ap = 4B, H, = 38.410 gauss 


cm 


For this, Py =0,112w. 


We shall now determine the resistance of the distributor 
coupling circuit by using the equation 

Ws 

Here, k,,=w,/w, is the ratio of the number of turns in the transmitting winding to the number of turns in the receiv- 
ing winding in the PD coupling circuits, Under optimum conditions with respect to power consumption, the coeffi- 
cient ky is assumed to be equal to 2. However, in order to reduce the so-called "reverse information,” it is con- 
venient to use larger k values, 


"Reverse information” can be characterized as the ratio of the transmission factors of the “reverse” and 
“forward” magnetic fluxes during the operation of a given magnetic element: 
AB 
where AB, and AB, are the magnetic induction variations in the subsequent and the preceding distributor elements, 
respectively, which are caused by the change in the operating element magnetic induction. 


We shall designate the reverse information coefficient by —. It can be shown that, in the absence of shunting 
rectifiers in the coupling circuits, this coefficient is determined by the value of kw: 


Ru = 


— 
—— 
— 
| — 
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(11) 


Thus, by choosing, for instance, k% =4, we can secure a sufficiently reliable PD operation (I= 1/16) without 
using additional shunting rectifiers, If k =4 and the number of turns in the receiving winding is = 50, we obtain 
Ry = 21.5 ohms from expression (12). 


The obtained resistance Ry must be several times larger than the forward resistance of the coupling circuit rec- 
tifier. For D-7 junction-type diodes, the forward resistance tf does not exceed 1 ohm under operating conditions, 
From this point of view, the obtained resistance R can be considered as acceptable. In the opposite case, different 
w, and ky values must be assigned. 


If, for any reason, the reverse information coefficient is not satisfactory, shunting rectifiers are installed, in 
which case the coefficient € is calculated from the equation 


By substituting Ry = 21.5 ohms, r¢= 1 ohm, and ky =4, we obtain k = 1/300, 


In order to determine the pulse source current and voltage, we shall compose an equivalent circuit for the 
operating PD cell (Fig. 8). The equivalent circuit is constructed in correspondence with the wiring diagram of a 
double-throw PD with shunting rectifiers for the limitation of “reverse information.” The forward resistances of all 
rectifiers are assumed to be equal to zero, In this circuit, Rc, is the resistance of the element magnetization reversal 
in operation: 


2B, P 
Wy 


Reg is the magnetization reversal resistance of the subsequent element: 


(13) 


R Pu 
HTN (14) 
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The power dissipated in resistor Rcg is equal to Py. 
In the equivalent circuit of Fig, 8, all resistances are reduced to the working winding, and, therefore, 
wy 1 
Roy= R 
w 


ww 


Ris = =) is the reduced resistance of the subsequent 
1 


element coupling circuit, By taking into account expressions (13 
— and (13), we obtain 
ly Ry 
k, —1 
6 ras The power dissipated in resistor Ry is calculated according 
4 * to the equation = (ky Rag = Rey (kw I) is the re- 
Fig. 8. duced resistance of the preceding element coupling circuit, The 
power dissipated in this resistor is calculated by using the equa- 


Ww 


tion is the reduced load resist- 


By using the above equations, we shall calculate the equivalent circuit parameters, For the determination of 
Rey, the number of turns wy must be given, In this, it is necessary to take into account the presence of the winding 
space as well as the fact that the required current from the pulse source decreases with an increase in wy. If w,=100, 
Roi = 28.7 ohms, 


ms, 
= 
nof 
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t= 
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Correspondingly, we determine the other parameters: 
Ri, = 1.8 ohm. NI = 5.40hm, Ni = 86. ohm. Nl. 1. 25 ohm. 


Since, with the exception of Ry, all resistances in the equivalent circuit as well as the power dissipated in these 
resistances are known, the current I through resistor Ry can be readily calculated, This current is necessary for the 
operation of the distributor cell in the circuit of Fig. 7. The current values in the equivalent circuit networks, re- 
duced to the primary * winding, are given below: 


i= 1A emp. = 0.0625amp, 


112 = = 0. 25amp, I, V5 = 0.0208amp. 


In order to check the calculations for the equivalent circuit, voltage Ua must be determined with respect to 


the known resistances and currents in the branches; this voltage must be of course equal for all the branches, In the 
case under consideration, Uap ~ 1.7 v. 


The total current is Ay = + t+ + 4, = 1-73amp. 
Assuming that the working pulse source voltage is equal to 24 v. we find the resistance of the distributor mo- 


tion circuit: 
Ry= = 13,20hm, 


In the latter equation, Ry = 0.645 ohm is the resistance of the equivalent circuit parallel branches, The resist- 


ance Ry comprises the ohmic resistances of all the Ww windings that are connected in series and the additional re- 
sistance Raq. 


Smsec 


Fig. 9. Pulse source, The core is made of Steel Fig. 10, Output pulse in the double-throw PD cir- 
10; do=60mm, din =40mm; W, = 1200 turns of cuit of Fig. 3 (the current through the Ry =40ohm 
0. 35- mm wire; C=6 uf. RA = 7 ohms. load). 


— 
Fig. 11. Oscillograins ot pulses in the single · throw PD Fig. 12, Output pulse in the single-throw PD cir- 
circuit (FigA), a) Voltage across capacitor Coon; b) cuit of Fig. 4 (voltage across the on capaci- 


voltage across the Rl 2300 ohms load; o) voltage across 
the Ry =6000 ohms load and across the CIA: capacitor, 
which is connected in parallel, 


tor for Ry =). 


— 


In conclusion, we shall give here the characteristics of some pulse distributors based on the circuits shown in 
Figs.3 and 4 and using E310 cores (Table 2), 


In all cases, 30x 20 * 6 mm cores are used, The rectifiers in the circuits are of the D7D type. 


The pulse distributors whose parameters are given in Table 2 allow considerable voltage fluctuations in the sup- 
ply sources (of the order of 50% of the rated value). 


The temperature tests, which were performed on a distributor based on the data in item 1 of Table 2,proved 
that the operation of this distributor is stable at temperatures attaining 100°C, 


As was shown by suitable experiments, the use of pulse distributors with magnetic cores made of Steel £310 in 
remote control devices with time separation of channels in using a synchronous 50-cps line as the supply source 
makes it possible to dispense entirely with intermediate amplifiers between the distributor and the output final units 


(MKU-48 relays, commutator incandescent lamps, etc.), which lowers the costs and makes the entire remote control 
device more simple. 
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ON THE DESIGN OF PUSH-PULL MAGNETIC AMPLIFIERS 


V. L. Benin 


(Khar*kov) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 3, 
pp. 393-399, March, 1961 

Original article submitted July 22, 1960 


The method of designing push-pull magnetic amplifiers presented in [1] is further developed to per- 
mit direct determination of the optimum parameters of these amplifiers without the need for succes- 
sive approximations, The accuracy is determined by the accuracy of the theory of linearized mag- 
netic amplifiers, 


In recent years, many methods of designing different types of magnetic amplifiers have been developed. 


The method given in [I] is one of the most rational engineering methods of designing magnetic amplifiers. 
The authors base this method on the theory of linearized magnetic amplifiers, and use the method of successive ap- 
proximations to design push-pull magnetic amplifiers with a complex load. 


In the first approximation, the optimum specific volume of a push-pull magnetic amplifier is determined from 
the minimum specific volume curves of inductor magnetic amplifiers, after which the magnetic amplifier param- 
eters are obtained with greater precision, according to [I]. 


The error of the first ar proximation can be significant, especially if the load is active and the current multi- 
ples are k <5, and thus one or two approximations must be made in order to get the required design accuracy, In 


| 
— — 
— 8 
— 
— 


this work, the design method given in [1] is further developed to permit direct determination of the optimum param- 
eters of push-pull magnetic amplifiers without successive approximations, The accuracy of this method is determined 
by the accuracy of the theory of linearized magnetic amplifiers. 


1. Specific Volume of the Steel in a Push- Pull Magnetic Amplifier 
The circuit diagram of a push-pull bridge magnetic amplifier is shown in Fig. 1. 


According to [1], the volume of the steel in one core of a push-pull magnetic amplifier is found from the for- 
mula 


10] sin 9%. +V sine ()- a 


H (By — By.) 
= 21 


H (By — BY) 


(2) 


The following designations are introduced in these expressions: P the apparent load power in volt-amps, Vcj~ 
the volume of one core of an inductor magnetic amplifier, in cm*, vi- the specific volume of one core of an induc- 
tor magnetic amplifier in cm*/volt-amp, Zir- load impedance in ohms, Ik-effective value of the current in the 
magnetizing inductor winding of a push-pull magnetic amplifier, in amps, Bx and B; the magnitudes of the alternat- 
ing flux density components in the magnetizing and demagnetizing cores of the amplifier, in gauss, H- the alternat- 
ing component of the magnetic field strength in the magnetizing core, in amps/cem, k=I/Ig=Hy/Hg~ ratio of the 
alternating operating currents in the magnetizing and demagnetizing windings of the cores (the current change multi- 


ple), Y—the phase angle between the voltage UH and current Ii in the amplifier load, and w= 2rf the angular fre- 
quency of the supply voltage. 


To determine the current Ix from the load parameters given in the design 


(Uy and I). we use the vector diagram of a push-pull magnetic amplifier shown 
in Fig, 2. 


When the resistances of the operating windings and core losses can be neg- 
lected, as is usually assumed in magnetic amplifier design, the phase angle be- 
tween the corresponding currents and voltages in the amplifier working windings 
(Ux, IK. Us, Ig) is equal to 90°. In view of this, the phase angle between currents 
IK and Ip is equal to the phase angle 8 between voltages Uy and U. 


It is easy to obtain the following equations from the vector diagram of Fig. 2: 
U2 = U*+ U2 20 U, cosa, (3) 

Un = U*? cosa, (4) 

4U* = Uz + 20, U, cos B, (5) 

= + 2. I, cos B: (6) 


The voltages across the operating windings of the magnetic amplifier can be expressed in terms of the core 
volume, the magnetic circuit parameters By, By, Hy, and the current Ix . Indeed, keeping in mind that 


Fig. 1. 


where w,, is the number of turns in one coil of the operating winding, and I and q are the mean length and cross 
section of one core, we find that 


* — ql, 


Vv 
U. = = @ByHy 
K 


Ä! 
NDS 
— 
j ' 


V 
10 = 7210". 
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Solving equations (3)-(8) simultaneously, and keeping in mind that 


(9) 
we get the equation 


1 ‘i (10) 
— B. Hu 10°" 2 = 0, 


from which we determine the current in the operating winding of the magnetizing core, Ix: 


(11) 
V By 
— 24 wal! +(e 


A = ＋ A, B. Hic 10". 


Substituting (11) into (1), we find the volume of the steel in one core of the push - pull magnetic amplifier, 
after certain transformations, 


1 


The specific volume of the steel in one core of the push-pull magnetic amplifier is found to be 
Ve 


20, 
Zulu 


Taking (2) into account, the specific volume of the steel in the core equals 


(14) 


B 5 
— K 
By, (14a) 


if the amplifier load is pure ly active (cos¢ 1). 

For a pure ly inductive amplifier load (cos Y = 0). 
251 
(14b) 
The lest expression fully agrees with the formula derived in [I]. 


The specific volume of the steel in the core of a push-pull magnetic amplifier is found with the aid of formula 
(14) from *he magnetization curves of a ferromagnetic material simultaneously by direct and alternating currents, and 
from the core steel specific volume curves of inductor magnetic amplifiers, In this case, as in [I], a current change 
multiple k must be chosen for each de value in such a way that the specific volume of the steel will be minimized, 
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(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(14a) 


2. Choosing the Optimum Alternating Current Change Multiple 


We use (14) to find the optimum ac current change multiple k. from the viewpoint of 
minimum amplifier weight, having approximated the magnetization curve by an analytic 
function, It is known that the magnetization curves of transformer steels and a number of 
other magnetic materials can be adequately approximated by the hyperbolic sine 


_ where-h is the instantaneous magnetic field strength in amps/cm, b—the instantaneous flux 


density in gauss, and q and 6 are coefficients characterizing the properties of the magnetic 
material, 


The quantities q and 8, computed from the magnetization curves of different materials 
given in [I]. have the following values: 


Material a, amp/cm B ,1/gauss 


£-42 0.12 4.2 
£-320 0.055 2.9+1074 
50-NP 0.005 5-10°¢ 


Let us determine the optimum current change multiple E as a function of the de magnetic field strength Hy 
for two extreme cases, cos =1, and cos = O. To do this, we will rewrite expressions (14a) and (14b) as dimension- 
less quantities, taking equation (2) into account also: for cos? =1 


le == 1— ’ Ma) 


for cos ¢ = 0 


BB,., 05 = BB,, 


The optimum multiple k can be found if n, Ok. and y are known as functions of the term ag. These dependen- 
cies are easily obtained from a normalized analytic expression of the simultaneous de and ac magnetization curves 
of the core with free, even, current harmonics, the fee curve be ing approximated by (15) [2]: 


ay = [1+ Ju. (20% — 1]. (21) 
71009 
where I (OK) is a Bessel function of the first kind, zero order, with a purely imaginary argument. 
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When aq >50, which la the case for transformer steels where Hq >5 amps/cm, we find, from (21), that 


Setting ag =0 in (21), we find 


H 
aK I. (23) 


and we get the equation 


“4 100 


from (22) and (23), from which the dependence aq/k= f(y) can be constructed for different values of OK. The fune- 
tions n= f (Qx), and ag/k=f (y), derived from (22) and (24), are plotted in Figs. 3 and 4. 


The optimum k for different ag is determined in the following manner, For 
207 a : each ay, and different k and OK. n and y are found from the curves of Figs. 3 and 4, 


and vz is found from expressions (16). The optimum ko, no, Ok e, and 50 are the val- 
— ues for which vj is minimum. 


lA Figure 5 shows the calculated dependencies of the optimum normalized param- 


eters of push-pull amplifiers as a function of the relative de magnetic field strength 
ag. Since the values of kp almost agree for active and inductive loads, the function 
Ro=f(a,) is shown only for cosy =1 in Fig. 5. It can be seen from the curves of Fig.§ 
that the optimum current change multiple k varies insignificantly over a wide range 


as of de magnetic field strengths, and when ay > 200, K% 5 (for £-42 steel, Hg>20 
‘ amps / em: for F-320 steel, Hd 11 amps / em: for NP-50 permalloy, Hq > 1 amp/cm), 
Fig. 3 When 50 < ay 200, the mean value of the optimum multiple can be taken as kg=4, 


The optimum ng over the whole range of ag changes from 1.4 to 1.65, which 
agrees with the conclusions in [I]. 


70 Since the optimum ko 5, the determination of the specific volume of the cor 
steel in a push-pull magnetic amplifier from the relations in CI] is in error by a sig- 


* 8 7475 nificant amount (more than 20%) when cos v =1; this error is eliminated by succes- 


sive approximations, 


It therefore seems desirable to do the calculations with the formulas given in 


Be | | this article, setting kg=4— 5. 


The optimum k for actual magnetization curves may differ somewhat from the 
/ values given in Fig. 5, but the specific volume of steel computed from (14) will dif- 
i 


fer insignificantly from the minimum. 


AAS 3. Concerning the Design of a Push-Pull Magnetic Amplifier 
A 


The design of a push-pull magnetic amplifier can be done analytically with 
7 satisfactory accuracy. Figure 6 shows the curve of the minimum normalized specific 
volume of the steel in the core as a function of the de magnetic field strength, T 
Fig. 4. curve can be plotted in absolute magnitudes by using known values of q and 6 (for 
example, see [3]), and all of the parameters found by the method given in [1]. One 
need only remember that in finding the number of turns in the operating winding, the current in it. Ix. should be de- 
termined from formula (11) of this article, and the supply voltage U in (3) and (4) from the formula 


＋ U2 


where Us and Ux are found from (7) and (8). 


2000 — 1. (29) 
— 
— 
— 
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Fig. 5. — cos = 1. cos¢ =0, Fig. 6. 


A bridge magnetic amplifier using Sh-32 cores made of £-42 transformer steel (a = 0,12amps/em, 6 = 
= 4,2° 10˙ gauss”) was designed analytically from the given dependencies; the active load was 200 watts for Hy= 
= 25 amps/cm, k=5, The difference between the computed and experimental values of load current was about ah, 
for nominal values of U and Hy. 


CONCLUSIONS 


1, A push-pull magnetic amplifier can be designed without successive approximations, and with an accuracy 
determined by that of the theory of linearized magnetic amplifiers. 


2. The optimum current multiple k=4— 5 for push-pull magnetic amplifiers, from the viewpoint of minimum 
weight, It is therefore desirable to design them from the dependencies given in this work, and not from the relations 
given in [I]. which require successive approximations for such values of R. 
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NOMOGRAM FOR DETERMINING PHASE BY MEANS 
OF AN ASYMPTOTIC LOGARITHMIC AMPLITUDE- 
FREQUENCY CHARACTERISTIC 


B. M. Menskii 


(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 3, 


pp. 400-402, March, 1961 
Original article submitted October 5, 1960 


It is known that a unique relation exists between the amplitude-frequency and phase-frequency characteristics 
in minimum phase systems and elements, The phase ¢ at frequency wecan be computed from the formula 


@ (00) = me 


where u = 1 r. di / du is the slope of the logarithmic amplitude · frequency characteristic, and c is a constant de- 
termined by the dimensionality of dA, /du and v. 


However, it is difficult to compute phase from equation (1), even with the use of an asymptotic logarithmic 
amplitude · frequency characteristic. Therefore, different methods of simplifying the calculation have been developed, 
Using the graph Incth|u/2|, which is given in [I]. it is possible to evaluate the phase approximately at chosen points 
[2]. Approximate formulas are used in [3] and in other works, A method of plotting the phase-frequency character- 
istic by using auxiliary graphs was developed in [4], The tables given in [5] can be used instead of these graphs. 
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Fig. 1. 
The phase at a chosen frequency can also be determined from an asymptotic amplitude - frequency character- 


istic by means of the nomogram given in Fig. 1. 
The nomogram is a graph of an auxiliary function 
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To find the phase at the chosen frequency w», one must place the nomogram on the asymptotic characteristic 


so that the arrow on the nomogram coincides with this frequency (Fig. 2). The phase in degrees is found from the for- 


value of S over the frequency range corresponding to the i-th asymptote; and n the number of asymptotes in the fre- 
quency range covered by the nomogram. 


plotted on a logarithmic scale. The nomogram is symmetric about the line w/w»=1, Computed values of 8 for 


where vi is the slope of the i-th asymptote in db/decade (taking sign into account); 


are given in the table, which was formed with the aid of the tables given in IJ. 


Fig, 2. 
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If the frequency w divides any asymptote into two parts, each of them must be considered as a separate asymp- 
tote, 


The phase can be found approximately by using a nomogram plotted over a range of two decades (as is shown 
in Fig. 1). To increase accuracy, one must use a nomogram plotted, by means of the data in the table, over a four 
decade range. 


Having found the phase for several frequencies, one can plot an approximate phase · frequency characteristic. 
Example. Let us find the phase at w»=8 from the asymptotic amplitude · frequency characteristic shown in 
Fig. 2: 
( = [0 (2.25 — 1.44)] + [— 40 (1.44 — 0)] + 
+ [— 20(4.32 — 0)] + [20(2.25 — 1.32)] = — 65.45. 


To check this value, we find that 
 (@o) = — 65,8° 


for the same frequency, by using the tables given in [5]. 
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ON THE METHOD OF GENERATING RANDOM PROCESSES 
WITH A GIVEN MATRIX OF SPECTRAL DENSITIES 
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A method of generating stationary random processes is generalized for cases when the spectral densities 
of the desired processes are not fractionally-rational functions of the quantity s= jw. The method is il- 
lustrated by examples of generating discrete random processes and continuous random processes with de- 
lay. 


Statement of the Problem 


In [I]. a method is described for generating continuous stationary random processes whose direct and reciprocal 
spectral densities are fractionally-rational functions of s= jw. The transfer functions of the generator filters (the so- 
called shaping filters) in this case are fractionally-rational functions of a complex variable. 


A method is described in this article for generating stationary random processes with any type of spectral density, The 
transfer functions of the shaping filters are determined in the same manner as in II]: however, they are not fractionally-ra- 
tionai functions of s in the general case. 


The generation of random processes by the method described in [I] is based on theorems about the properties of 
spectral density matrices (see the appendix in [1]), These theorems, evidently, are true for spectral densities of any 


arbitrary type in the class of stationary random processes, since their proofs are based only on the properties of Her- 
mitic matrices, 


The formulas defining the transfer functions of the shaping filters also hold for cases other than when they are 
fractionally-rational functions of 3. The restriction in [1] was made only so that the transfer functions of the shaping 
filters would be fractionally-rational functions of the complex variable 3. 


By discarding this limitation, one can use the formulas given in [I] to find the transfer functions of shaping fil- 
ters for generating random processes with arbitrary spectral densities, 


Continuous Random Processes 


All the statements and formulas of reference [1] remain true for any kind of spectral densities, Therefore, they 
are not stated here, but only examples of their application are given. 


Example 1, Suppose it is necessary to generate continuous random processes ui and ug, whose spectral densities 
have the form 


i—s* 10 — s* i+s 
Vn (e) = » (s ) — 1G (8) 2e 


From formulas (11) and (12) in [II]. we get 


1 1 
Yn (s) = 2 1 Hi (0). Yu (s) = e (i 8) Ay, (s). 
The function H,(s) must be defined in such a way that Tu and Vn will be transfer functions of stable filters, It 
follows from this condition.that 
Hi (s) = ge ° 
By substitution, we get 


rug Yu) 


Now we will find the auxiliary function Ge: 


(8) = (8) — Vu (8) Yin (— 8) = 
Whence 


3 
Yu) = +38) 
A block diagram of a generator which generates the desired processes uz and ug is shown in Fig. 1. 


Random Sequences 


The spectral densities of random sequences (discrete processes) are ordinarily expressed as functions of the 


quantity z=e* (s= jw), which corresponds to relative time (the pulse repetition period is T= 1). 


If 01 (2) is the direct spectral density, a function Aj(z) always exists such that 


§G,, (2) = A; (2) A, (-). 
This equality is analogous to formulas (1) and (16) in [1], since z~*= e*§, 
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The generator diagram is completely analogous to the diagram in Fig. 1 of reference [I]. where Qi are genera- 
tors of non-correlated “white” random sequences (discrete processes), and Fix are discrete filters. 


2(t-s) 


rig. 1. rig. 2. 


Ultimately, one can arrive at the conclusion given in [1]. In this case, the transfer functions II (z) of the dis- 
crete shaping filters are found from the formula 


( H. (2). 


The auxiliary spectral densities Gin are found from formula (13) in [1], and Hi (2) are transfer functions of dis- 
crete phase filters; in general form 


Hy) = 725. 


If the spectral densities of the required random sequences are fractionally-rational functions of z, the functions 
Gin: Aj, Hj, and consequently ix - will also be fractionally-rational functions of z. The transfer functions of the 
phase filters in this case assume the form 


— (m>0, 11.2. ). 


12 — 
Example 2. Let it be necessary to generate discrete random processes ui and uz, whose spectral densities are 


» — 1.5 + 3.252 — 4.52? 
(2) = 1.252 + 0.52 
—0.8 + 2.22 — 1.52 
(2) = + 1.362 — 0.622" 
—0.4 + 2.82 — 0.42" 
"Gas (2) = + 1.362 — 0.62 


One can obtain 


— — 0.81 — 1.52 


as transfer functions. 
We find, from the condition that In be a transfer function of a stable filter, that, 


Ay (s) = 


z—1.5 z— 0.8 
vn (2) = 70405)’ rn (2) = — 0.60 


The auxiliary function C (2) has the form 


+ 1.162 4s? 
= 


ͤ— — 


so that 


(2) = 56. 005 


A block diagram of a generator which generates the required processes ui and uz is shown in Fig. 2. 


CONCLUSION 


A generator of any number of stationary random processes whose spectral densities are fractionally-rational 
functions of s= jw was described in [1]. This article gives a generalization of this method for cases when the spectral 
densities of the required processes have arbitrary form. In particular, it is possible to generate discrete random proc- 
esses (random sequences) by means of this method, by using digital computers, for example. 
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THE POSITIONING OF SENSING DEVICES IN A CONVERTER 
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Different methods of positioning sensing devices which operate on the V-reading principle in con- 
verters of angular rotation or linear displacement into digital form are examined, The use of these 

methods eases the construction of the converter, A method of deciphering the codes read off is de- 
scribed, 


False codes can arise in converters of angular rotation or linear displacement, operating on the reading princi- 
ple, into digital form. This occurs at scale points where the code changes by more than one digit, because of the 
finite size of the sensing devices, and the variation in their physical characteristics [I]. 


Certain methods of eliminating false codes exist, Converters are examined in this article which operate on the 
principle of a logical choice of the sensing devices, or according to the V-reading principle as it is frequently called. 


According to this method, two sensors, placed symmetrically about a reference line, are put in each digit posi- 
tion, except zero; their orientation reminds one of the letter "V" (Fig. 1), The zero digit sensor Vy is placed on the 
reference line, The displacement of the sensor relative to the reference line will be called Aj, and the size of the 
i-th digit, xi. The Vi sensors constitute the left group, and the Vi- the right group. The striated portions of the scale 
designate code "0", and the clear portions - 1“. Scale movement along the o arrow increases the code read off, 


The choice of sensors, in reading the code, is made according to the following logical rule. If the zero sensor 
Vo reads "0", we choose the sensor in the right group, Vj, for the first digit xi. If the zero sensor Vy reads 1.“ we 
choose the sensor in the left group, Vi, for the first digit x1. The figure read by the first digit sensor chosen is the 
criterion for an analogous choice of the sensor of the second digit x2, etc, 


The decoding operation, in connection with this type of logic, is done in a digital computer (Dig. Comp.) ac- 
cording to the diagram shown in Fig. 2. This circuit consists of two input registers (right and left), and two logic 
cells Inh. and "and." The registers consist of logic cells, each of which delays the information fed in for one operat - 


2 
| 
⁊ 
| 


ing cycle of the Dig. Comp. The logic cell Inh. inhibits when information is fed to the inhibit input. The "and" 
gate has an output only when both inputs are excited simultaneously. 


The decoder works in the following way, The information in the converter is read out on command simultane- 
ously at all digits, and is written in the appropriate cells of the right and left registers. Then this information is 
shifted from one cell of the register to the next after each operating cycle of the Dig. Comp. In this manner, the in- 
formation from the converter is fed in the form of a sequential code to the logic cells. These cells do the decoding, 
This kind of decoder makes it possible to decrease the required number of logic cells, and to decrease the number of 
wires needed for transmitting information to the Dig. Comp, Let us examine the effect of the displacement of the 
sensing devices relative to the reference line on the performance of the converter. Converters having their sensors dis- 
placed by A; =0.25 nN. where n= 1, 2, 3. , are inoperative, since the code changes simultaneously by more than 
one digit at certain scale points, Also, converters will not operate reliably if their sensors are placed within a certain 
zone near these points, These dead“ zones are approximately two or three times wider than the sensitive zone of the 
sensor, The remaining zones are effective, If the sensors are placed within the effective zones, they will be op 
tive. The location of the dead“ and effective zones is shown in Fig. 3. 
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The center lines of the effective zones are the opti- 


mum locations for the sensors; their coordinates are equal 
to 0,125xj, 0.375x;, 0.625x;j, 0.875x;, etc. The effective 
zone centered around the 0,125xj; line will be called zone I. 
zone 2 has the 0,375x; center line, zone 3—0,625x;, and 
zone 4—0,875x;. 


Figure 1 shows the locations of the sensors in zone 1 (Aj=0,125x;). The direct code is used for this case (all 
digits start from the dark areas), and decoding is done according to the diagram shown in Fig. 2. 


The location of the sensors in zone 2 is shown in Fig. 4. An inverse code is used in this case for correct opera~ 
tion of the converter, The inverse code scale differs from the direct code scale in that white areas are replaced by 
dark, and dark by white, Decoding is done according to the diagram shown in Fig. 5. 


For direct connection of the sensors, the left sensor Vi is connected to the left register, and the right sensor Vj 
is connected to the right register, For inverse connection, the left sensor is connected to the right register, and the 
right sensor to the left register. The zero sensor is connected to the right register in all cases, 
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The location of sensors in zone 3 is shown in Fig. 6. The inverse code is used in this case. Decoding is done | 
according to the diagram of Fig, 2, Direct connection of the sensors is required in this case. | 


Figure 7 shows the location of sensors in zone 4. The direct code is used, and decoding is done as in Fig. 5, In- 
verse connection of the sensors is required here, 


The table shows the varieties of code, and the method of connecting the sensors, depending on their location, 


Zone number Relative displacement Code Sensor connection | 
1 : 0.125 direct direct 
2 0.375 inverse inverse . 
3 0.625 inverse direct 
4 0.875 direct inverse 


One should keep in mind that in V-reading type converters sensors for each individual digit can be located in 
any zone which is most convenient for structural reasons, The code and method of connecting the sensors for each 
digit depend on which zone the sensors are in, in accordance with the table, 


AVA 


Witt: Wh —— 


rig. 7. 


CONCLUSIONS 


The sensing devices in a V-reading converter of angular rotation of linear displacement into digital form can 
be placed at any part of the scale, except dead“ zones, by placing two sensors symmetrically about a reference line 
for each digit in the code. The zero digit sensor is an exception, There is only one sensing device, and it is placed 
on the reference line, 


For each digit except zero, there are four methods (or zones) of placing the sensors, depending on the amount 
by which they are displaced from the reference line. The coding and method of connecting the sensors to the regis- 
ters depend on the zone of location. 


Construction of a converter is a complex problem if the sensors are placed only in zone 1, The use of all zones 
for locating sensors significantly eases this problem, 
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THE CALCULATION OF A TWO-TERMINAL NETWORK 
WITH ONE SEMICONDUCTOR THERMISTOR OPERATING 
UNDER LINEAR CONDITIONS 


N. P. Udalov 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 3, 
pp. 409-412, March, 1961 

Original article submitted May 23, 1960 


The article presents the calculations for a two-terminal network with a semiconductor thermistor, 
operating under linear conditions, which can be applied in circuits with temperature compensation 
and microthermostats, 


The calculation of a circuit for temperature compensation of a linear resistor which has a positive temperature 
coefficient can be reduced to the determination of the parameters of the elements in a two-terminal network which 
has a temperature characteristic approaching linearity, and an increase in resistance, during ambient temperature 
changes within given limits, equal to and opposite in sign from the corresponding increase in the resistance of the com- 
pensating element. 


As is known, these requirements are satisfied by the two-terminal network shown in 


I, Rk Fig. 1, where R is a semiconductor thermistor (SC T), 11 is a resistor which changes very 
ee little with temperature and is considered of fixed value in the calculations. The tempera- 
| re | ture characteristic of such a two-terminal network in its best form approaches linearity if, 
h for temperatures T, = Tmin + O. IA T and Tz = Tmax — C. IAT. where AT = Tmax Tmin, 


Fig, 1. Two-terminal it intersects with a given linear characteristic, and for Tay = A157 u has a point of 


network with one fixed discontinuity [1], The latter condition is met when the following equality is observed: 
2 B — 27. 
— 


"= B+ Pav’ (1) 


where Ray is the SCT resistance at a temperature Tay. B is the coefficient of the SCT temperature characteristic, 
the value of which is determined from the values of the SCT resistances Ri and N at temperatures TI and Tz: 


T,T R 
B= In (2) 


Inasmuch as the SCT operates under linear conditions in the circuits under examination, I. e., with currents 
which do not produce noticeable heating of the active part of the SCT, the quantities Ri. Ray, and N may be ex- 


pressed by the resistors of the two-terminal network, rk:, rx. ay» 20d tke, Corresponding to the temperatures TI. Tay. 
and T;: 


Fav = ay’ (3) 


After substituting expression (3) in (1) and (2), we obtain 


71 
2Tay 71 
B T,T: — 
(5) 


— 
— 
— 


If the linear temperature characteristic of the two-terminal network is given (Fig. 2,a), then in expressions (4) 
and (5) there remain only two unknown quantities B and ri, which can be determined simply by constructing graphs of 
the functions (4) and (5) with values for r1 in the interval 


Pig k. av- (6) 
The coordinates of the intersection points of these graphs (Fig. 3) correspond to the values B and ri, which we are 
seeking. 
The temperature characteristic of the SCT is defined by the 
r, ohms expression 


The SCT parameter Ray inserted here may be found from 
one of the expressions in (3), 


f The temperature characteristic of the compensator (Fig. 2. U 
* may be constructed by substituting in the expression 
mR 


the values of R found from the SCT temperature characteristic for 
a series of temperature values, 


Example 1, To calculate the compensator for a copper re- 
sistor gy =300ohms (for T,y=273°K) if Tmin=223°K, Tmax 
= 323°K, 


We find AT = 100°, T,=233°K and the corresponding copper 


| 
| resistor tei 249.5 ohms, T2 313 K and tc, = 350.5 chms. The in- 
f | as crease in the copper resistor Ar= 101 ohms, Let us assume that the 
mi ae RWS linear temperature characteristic of the compensator satisfies the 
| 4 1 condition |r — cuil Ar (Fig. 2. a), in accordance with which ru 
a) | = 121 ohms, ay = 70.5 ohms, rpg 20 ohms. 
2 Tua 063, 5803 K Substituting the known values in (4) and (5), we obtain 
Fig. 2. Temperature characteristics of a to- te- B = 546 wo B = 911.48 in 
minal network: a,b apply to example 1;c, d- . 
to example 2. Solving this system of equations graphically (Fig. 3), we find 
(2) t,= 125.6 ohms, B= 4500°K (the values for ti were taken in 
the interval 121 <r, < 141 ohms), The SCT temperature character- 
aw | istic in conformity with (7) is written in the form 
ye 4500 4500 
R = 160 exp (F — 273 
22. The temperature characteristic of the compensator is con- 
(3) \ "J structed in Fig. 2. b. 
— 
— In the example under consideration the line ar characteristic 
of the compensator was given in accordance with the requirements 
t for minimum power consumption in the compensator (ideal from 
(4) on ae are this viewpoint would be a compensator with zero resistance at 


Tmax) 
Fig. 3, Graphical solution of a system of equa- 
(5) tions (applying to example 1). Let us assume that the given characteristic corresponds to 


Fig. 2. c. In this case, with the same conditions of the problem, 


| 
| 
| 
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we shall obtain; r,; 496 ohms, B=1100°K, Ray=14000hms, The temperature characteristic of the compensator is 
shown in Fig. 2,d. Comparing curves b and d in Fig. 2, it is easy to note that the gain in power results in a decrease in 
quality of compensation, which may be determined by the maximum value of the relationship A/ te, ay (here A is 
the deviation of the compensator characteristic from the given linear relationship, andr, ay is the resistance of the 
compensating element at a temperature T,,). This value is 3.3% in the first variant of our problem, and 0. 5 / in the 
second, 


It appears possible to find a quantitative relationship between the quantity A/ fe. gy and the relationship of 
k. ay to fc, ay which permits the use of a linear characteristic for the compensator based on the given accuracy of 
compensation, 


In practice, it is not always possible to select SCT with parameters coinciding with the calculated values, In 
that case, any type of SCT may be selected which satisfies the condition |R,—R,| > Ar. The compensating circuit will 
then contain two fixed resistors r, and 12 (Fig. 4), and the calculation problem is reduced to the determination of these 
resistances, For its solution we shall make use of the condition specifying the existence of a point of discontinuity on 
the temperature characteristic of such a two-terminal network [1]: 


B — av 
Rav (9) 


and the condition of equality of the increment in the compensator resistance of a given value rky tr = Ar. Substi- 
tuting here 
(Re + ra) 


(Ai + ra) ri 


and solving the latter equation with respect to r,, we have 


Ar + + Ra) (ri + ra) + + (10) 


Example 2, To calculate the compensator in accordance with conditions in example 1, when the given SCT has 
a temperature characteristic R= 1.57 exp(1890/T). 


We find Rgy=15700hms, From expression (9) ri rz 870 ohms, and from 
expression (10) 11 450 ohms, i. e., 12 420 ohms. The temperature characteristic 


I, 
— of the compensator is shown in Fig, 2. d. 


13 The design of a simple microthermostat intended for the thermostatic con- 

Ih trol of separate circuit elements [2] involves the determination of the SCT param- 
Fig. 4. Two-terminal network eters which, with a change in temperature of the medium within given limits, will 
with one SCT and two fixed re- produce a change in current through the resistance ri which conforms to a law ap- 
sistors. proaching linearity, The resistor ri (Fig. 1) in this case consists of the heater coil 
N located on the thermostat element. The resistance ri and the current through it, 
Ini and Ihz at temperatures TI and T, are given. The current Io, before it splits 
up, is held constant by means of any suitable regulator. 


711 


After constructing the volt - ampere curve of resistor ri, we shall mark points on it 
which correspond to the given values of Ip; and Ip, (Fig. 5). If we assign a value 
for the current Io, then from Fig. 5, the SCT parameters Ri and Ra will be deter- 
mined from the expressions 


R, =7— R= 


15 


It is quite obvious that the current Io may have any value in the interval 
Ih; < Iy< and should be taken as small as possible. Since I, is always greater 


Fig. 6. The, determination.of the chan Inn: then the value of theSCT resistance at a temperature Tz is also limited: 


SCT parameters in a microther- 
mostat circuit, Th, 
~ "he 


1111 
— 


& 


0) 


If the SCT is given, i. e., the values of its resistances at temperatures TI and Ta are known, then it may be used 
together with an additional resistance 12 (Fig. 4). The unknown quantities in this case will be the resistance ry andcur- 
rent Ig, To determine rz, we may make use of the condition of constant current in the common part of the circuit, 
which, for temperatures TI and Iz, will be 


Ini Ther 
1 1 =], 
m te 
from which 
2 ‘ Im — The 
The current in the common part of the circuit will be / 
I 
ha”? 
Ina + 


The above method of calculation is based on the assumption that the SCT operates under linear conditions, 
However, in microthermostat circuits where the current through the SCT varies within quite wide limits, considerable 
heating may occur in the active part of the SCT. Therefore, after determining the circuit parameters, this same prob- 
lem must then be solved by the graphical method [2]. 
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